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ABSTRACT 


•  • 


The  most  general  numerical  formulation  that  can  be  made  of  the  aero¬ 
dynamic  lifting  problem  is  a  matrix  of  aerodynamic  influence  coefficients 
(AlCs)  that  relate  th^  control'C collocation)  point  forces  (the  local  surface 
integral  of  the  pressure)  to  the  control  point  downwash  velocities  (or  the  control 

point  deflections,  who*8e  substantial  derivatives  are  the  downwash  velocities). 

•  » 

The  AlCs  are  necessary  in  any  aeroelastic  analysis  by  collocation  methods. 
Furthermore,  the  AlCs  offer  a  large  computational  advantage  in  aeroelastic 
analyses  by  modal  methods  since  the  AlCs  need  be  derived  only  once*and  the 
generalized  forces  follow  from  multiplications  with  the  modal  matrix;  any 

approach  that  derives  generalized  forces  directly  iV inefficient  to  the  extent 

*  •  •  • 

that  the  aerodynamic  analysis  must  be  repeated  for  every  revision  of  the 

modal  set.  *  **.^  •*  ,  * 

•  *  *  •  *  • 

•  •  •  • 

^  .Downwash  is  considered  to  arise  from  either  motion  or  a  gust,  and  the 
unsteadiness  is  considere'd  to  be  oscillatory,  transient,  or  static  (a  limit*of 

*  y  •  •• .  ,  ,  ^  • 

the  oscillatory«case)i  The  oscillatory  AlCs  are  defined  l5y  .*  *  * 
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for  motion,  and  by 
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•  • 
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for  a  gust.  The  static  AlCs  are  defined. by  °  “ 

o 

fsl  =<‘is/:c)[cjj.4jv 


The  use  of  AlCs  thus  defined  is  illustrated  in  collocation  formulations  of 
several  aeroelastic  problems,  and  the  use  of  the  AlCs  in  analyses  by  modal 
methods  is  discussed.  Several  possibilities  for  the  definition  of  transient, AlCs 
are  suggested  in  a  preliminary  review  of  the  transient  problem. 


•11 


•  % 


A _ m. 


The  status  of  AlCs  is  discussed  for  surfaces  of  various  planforms  by  two- 
dimensional  (i.e.  ,  strip)  theories  and  by  three-dimensional  theories,  for 
slender  bodies  and  for  panels,  throughout  the  range  of  Mach  number  of  current 
interest.  A  unified  lifting  surface  theory  approach  to  AlCs  is  demonstrated  in 
the  subsonic-transonic- supersonic  Mach  number  regimes  based  on  a  critical 

review  of  current  techniques.  A  systematic  approach  to  interference  effects  is 

* 

s^gested.  A  critical  assessn^ent  of  transonic  and  hypersonic  problems  is 
presented,  and  the  linear  approximations  of  slender-wing  theory  and  piston 
theory  are  discussed  in  detail  to  illustrate  the  AIC  development.  A  result  is 
derived  for  the  AlCs  from  slender  body  theory,  and  the  significance  of  fir  st¬ 
and  second-order  slender  body  theory  is  reviewed.  The  AlCs  for  a  two- 
dimensional  panel  are  given,  and  the  problem  of  adapting  three-dimensional* 
aerodynamic  theory  to  the  finite  panel  problem  is  mentioned.  *  The  oscillatory 
AlCs  for  a  hg-rmonic  gust  are  developed  from  incompressible  strip  theory,  and 

a  method  of  extension  of  the  AICs*for  motion  to  obtain  the  gust  AlCs  from 

•  •  ••  •  ,  •  * 

,  lifting,  surface  theory  is  shown.  A  discussion  is  given  of  current  investigations 

r  *  •  •  *  * 

of  possible  formulations  of  Jransient  AlCs  and  of  related  computational  prob- 

•  .lems.  An  empi/ical  weighting  matrix  that  permits  adjustment  of  the 

theoretical  AlCs  to ‘incorporate  Jthe  ^results  of^  experimental  measurements  is 

suggested,  and,  finally,  a  number  of\ec6mmendations  for  further  irivesti- 
•  •  •  *  •  * 
gations  and^experimeRtal  correlation  are’made.  .  t  •  ,  •  * 
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SYMBOLS 


nm 


,  b 

^  r  o 


Cp.  Cpi 


Element  of  matrix  relating  pressure  loading  function 
amplitude  to  downwash  at  point  i 

Element  of  flexibility  matrix;  speed  of  sound;  panel  dimension 
Amplitude  of  deflection  mode 

Amplitude  of  pressure  loading  function 

Elernent  of  integration  matrix  for  slender-wing  theory 

Element  of  integration  ma^rix  relating  control  point  forces 
to  pressure  loading  function  amplitudes 

Local  semichord;  constant  in  exponential  approximation  of 
aerodynamic  lag;  panel  dimension 

•  • 

Reference  semichord,  and  root  semichord,  respectively 


Skin  friction  coefficient  for  compressible  and.incompressible 
flows,  respectively 

Element  of  AIC  matrix  for  a  harmonic  gust 

Element  of  oscillatory  AIC  matrix  for  motion 

•  • 

Element  of  static  AIC  matrix 


•  o 

Pressure  coefficient 


Lift  coefficient;  is  lift  curve^slope 

o 

Root  chord;  panel  span;  c  is  control  surface  chord;  ~c  is  mean 

a 

aerodynamic  chord  • 


Element  of  matrix  relating  downwash  at  point  i  to  pressure 
loading  function  amplitude  *  * 


Dista*nce  between  forward  and  middle  control  points;  diameter 


rs 


F 

F 


vN^'  N 


vn 


g 

g^%) 


H 

H 


,H  .H 
mn  r  rs 


h 

I 


I  .  J 

n  n 


I*' 

r 

i 

K 


K 

K 


mn 

ij 

mn 


mn 


Amplitude  of  deflection  mode 
Control  point  force 

Downwash  factor  for  concentrated  chordwise  replacement 
load 

Element  of  matrix  relatinc  AP  to  a 

®  V  nm 

Structural  damping  coefficient 

^Spanwise  loading  function 

Boundary  layer  parameter,  H  =  6  /O 
Assumed  deflection  functions 

9 

Surface  deflection 

Element  of  unit  matrix;  momentum  of  cross-flow  virtual  mass 

•  •  • 

Slendei^body  spanw^ise  integrals;  piston  theory  thickness 
integrals  ^  • 

Element  of  deflecfeon  interpolation  matrix 

Imaginary  unit  •  •  * 

Kernel  of  downwash  integral  equation;  constant;  K0  is 

nonsingular  numerator  of  subsonic  kernel  function;  Kq  and 

•  •  ^  b  • 

K  are  hypersonic  similarity  parameters 

Slender-wing  theory  surface  integrals 

Elefnent  of  slender-wing  theory  aerodynamic  inertia  matrix  ^ 

Local  reduced  frequency,  k  =  b/V 
Reference  reduced  frequency 

Lift 

•  . 

Oscillatory  lift  coefficients  * 


viii 


Jt 


Surface  semispan  (in  Section  11.  B.  2,  otherwise  s  is  used) 
Chordwise  pressure  loading  functions  ♦ 

Ma'ch  number;  element  of  mass  matrix;  pitching  moment 
Oscillatory  pitching  moment  coefficients 

Element  of  slender-wing  theory  aerodynamic  spring  matrix 


Number  of  downwash  interference  fields 


p.Ap 


Element  of  slender-wing  theory  aerodynamic  damping  matrix 
Normal  vector 

■# 

Concentrated  chordwise  replacement  load 

Pressure 

Generalized  force 

Dynamic  pressure;  q  is  n*th*generalized  coordin<ii.c 
•  Reynolds  number 

Body  radius  in  cylindrical  coordinates  • 

Radial  cylindrical  coordinate  •  • 

Surface  reference  area;  body  cross-sectional  area 

•  •  • 

Incremental  are*a  surrounding  control  point 

•  •  ^  • 

Surface  semispan  (except  in  Section  II.  B.  2  where!  is  used) 

Hinge*moment;  representative  decay  time;  temperature, 

is  wall  temperature,  T^,.  is  mean  boundary  layer  temperature 

♦  *  *  * 
Oscillatory  Jiinge  moment  ct)efficients 


Time 


Streamwise  velocity  in  boundary  layer 
Fred  stream  velocity  * 


ix 


o 


Aq.  Aj 


Nonhomo  geneou  8  functions  in  second-order  slender -body 
solution 

Complex  eigenvalue,  X  ”  g^st  wave  length 

Coefficient  of  viscosity 


Dimensionless  coordinates 


Atmospheric  density 


*(k) 


Ratio  of  semispan  to  root  chord 

Dimensionless  time;  real  time  measured  in  inertial  reference 
frame 

Potential  functions 

Sears  incompressible  harmonic  gust  function  * 

Hypersonic  viscous  interaction  parameter;  pressure  potential 
function  for  flows  with  vorticity 

Second -order  potential  function 

Potential  function 
Frequency 


Subscripts  and  Supercripts 


Evaluation  at  body  surface;  refers  to  basic  weighting  factor 


,  Experimental 


Fuselage 

Evaluation  in  fluid 


Gust 


(  )‘ 


Index  of  downwash  control  point  grid 


Index  of  deflection-^oad  control  point  grid 


Refers  to  deflection  mode  of  amplitude  a 


xi 


(  ) 

<  >: 
(  ) 

(  ) 

(  ) 


nm 


s 

(t) 


(  ) 

’(  ) 

• 

(  ) 


W 


0 


1 


(  ). 


n 


(  ) 


Refers  to  pressure  loading  coefficient 
Rigid;  refers  to  redundant  weighting  factor 
Sta'tic;  shock  location 

Theoretical 

Partial  derivative^with  respect  to  x;  (  (  )^,  ( 

denote  differentiation  with  respect  to  y,  z,  and  t,  respectively 

Wing 

Free-stream  condition;  reference  value;  axial  flow  or 
symmetrical  potential  ^ 

Cross  flow  or  lifting  potential;  condition  upstream  of 
shockwave 

Condition  downstream  of  shockwave 
Complex  amplitude 
Derivative;  disturbance  value 


Matrix  Notation 


Square,  rectangular,  row,  or  diagonal  matrix  ^ 
Column  matrix 
Transposed  matrix 
Inverse  of  matrix 


I.  INTRODUCTION 


The  essence  of  the  aerodynamic  lifting  problem  is  the  relationship  between 
the  distribution  of  lifting  pressure  (the  pressure  differential  between  the  upper 
and  lower  surfaces)'  and  the  surface  distribution  of  downwash  velocity.  The 
mathematical  formulation  of  the  lifting  problem  leads  to  an  integral  equation  for 
the  lifting  pressure.  The  most  general  numerical  formulation  that  can  be  made 
?rom  the  aerodynamic  integral  equation  leads  to  a  matrix  of  aerodynamic  influ¬ 
ence  coefficients  (AICs)  that  relate  the  control  (collocation)  point  lifting  pres¬ 
sures  (or  the  control  point  forces*  which  are  the  local  surface  integrals  of  the 
pressures)  to  the  control  point  downwash  velocities  (or  the  control  point  de¬ 
flections,  whose  substantial  derivatives  are  the  downwash  velocities).  For 
aeroelastic  analysis,  AICs  that  lead  to  control  point  forces,  rather  than  pres¬ 
sures,  are  more  useful.  In  the  problems  arising  from  atmospheric  disturbance 
(e.  g.  ,  a  gust),  the  AICs  may  relate  the  forces  to  the  downwash  directly;  in  the 
problems  arising  from  motion  of  the  surface  (e.  g.  ,  lift  effe^ctivene ss,  flutter), 
the  AICs  should  relate  the  forces  to  the  deflections  or  the  streamwise  slopes.  In 

the  present  treatment,  we  prefer  the  deflections  rather  than  streamwise  slopes 

•  • 

as  the  generalized  coordinates,  because  the  deflections  have  a  more^general 
meaning  on  a  cambered  surface,  and  deflection  influence  coefficients  are  more 
readily  obtained  from  a  structural  analysis  than  slope  influence  coefficients. 
Accordingly,  we  have  defined^  a  matrix  of  static  AICs  fcfr  quasi- static  motion  by 


=  (qS/c)  [Cj|h 


We  have  defined  *  a  (complex)  matrix  of  oscillatory  AICs  for  harmonic 
motion  by  • 


f|  =pcoVs[cJ|h) 


We  have  defined"^  a  (complex)  matrix  of  oscillatory., AICs  for  a  harmonic  gust  by 


F  =  p  VW  b  s  C  w  /W 
g  ^  g  r  [  gj  g  g 


We  note  that  the  definitions  of  Eqs.  (l)to  (3)  are  completely  general,  being  equally 
applicable  to  two-dimensional  (i.  e.  ,  strip)  or  three-dimensional  aerodynamic 
theories.  We  also  note  that  the  forces  given  by  Eqs.  (1)  and  (2)  must  be  equal  in 
the  limit  of  zero  frequency 


(4) 


which  leads  to  the  relationship  between  the  static  AICs  and  the  oscillatory  AICs 
for  motion  *  • 


f  hs]  =  [s]  <5* 

Hence,  the  static  AICs  need  no  special  consideration  in  the  subsequent  discus¬ 
sion.  We  do  not  wi^h  to  minimize  the  importance  of  the  many  developments  in 
4- 12  * 

the  static  context.  However,  we  shall  concentrate  on  the  unsteady  formu¬ 

lations  because  of  the  recent  number  of  significant  advances  in  unsteady  aero¬ 
dynamic  theory,  while  recognizing  that  there  are^some  problem  areas,  as  we 
shall  point  out  in  our  review  of  piston  theory,  that  await  the  steady  solution  be¬ 
fore  the  imsteady  one  is  known. 

We  may  Semonstrate  the  general  usefulness  of  AICs  by  considering  the 
collocation  formulation  of  a  number  of  the  aeroelastic  problems.  For  the  illus¬ 
tration,  we  choose  the  problems  of  static  lift  distribution,  flutter,  and  the  fre- 
quency  response  analysis  required  for  structural  fatigue  analysis  under  random 
atmospheric  turbulence.  For  simplicity,  we  assume  no  rigid  body  degrees  of 
freedom.  Our  basic  equation  relates  the  structural  deformations  to  the  applied 
forces  through  the  structural  influence  coefficients  , 


For  the  static  aeroelastic  analysis,  the  external  forces  are  given  by  Eq.  (1) 
where  the  deflection  mode  is  found  from  the  sum  of  the  initial  deflection  mode  of 


the  surface  considered  as  rigid  and  the  deformation  mode 


(7) 


Substituting  Eq.  (6)  into  Eq.  (7)  and  the  result  into  Eq.  (1)  permits  solution  for 
the  external  loading: 


pal  =  [A!  Pr|  <8) 

where 

p^l  =  (qS/c)[Cj  |h^|  •  (9) 

and 

[A]  =  ^[I]  -  (qS/c)  [Cj^J[a]y*  (10) 


A 


From  Eq.  (8)  not  only  may  the  structural  loads  be  calculated  (i.  e.  $  sheart 
moment,  and  torque),  but  also  the  aeroelastic  stability  derivatives  may  be  eval¬ 
uated  (e.  g.  ,  the  lift  curve  slope  is  given  by  =  [l]  jF^j/qSor  when  we  take 

I  j  =  the  case  of  Asymmetrical  angle  of  attack). 

We  next  consider  the  flutter  problem.  A  gutter  analysis  is  usually  carried 
out  ^y  assuming  a  harmonic  solution  and  seeking  the  combination  of  speed  and 
frequency  for  which  harmonic  motion  actually  exists.  To  this  end,  an  artificial 
structural  damping  is  introduced  into  the  system  to  sustain  the  assumed  har¬ 
monic  motion,  and  the  flutter  point  is  found  when  the  actual  damping  and  the  re¬ 
quired  damping  are  equal.  The  total  applied  loads  require  the  addition  of  the 
inertial  loads  to  the  oscillatory  aerodynamic  loads  of  Eq.  (2),  which  for  the  as¬ 
sumed  harmonic  solution  leads  to 


3 


1 


1 


Combining  Eqs.  (6)  and  (lib),  and  introducing  the  artificial  damping  as  a  complex 
divisor  (a  complex  factor  to  the  modulus  of  elasticity)  to  the  flexibility  matrix, 
leads  to  the  formulation  of  the  flutter  eigenvalue  problem 

2 

I’^fl  =  |*'f|  • 

2 

The  eigenvalues,  \  =  +  iX.j  =  (1  +  ig)/w  i  of  Eq.  (12)  may  be  found  by  complex 

matrix  iteration  from  which  we  find  the  required  damping 

g  =  .  (13) 

the  frequency 

m  =  1/YX^  .  (14) 

and,  since  the  formulation  of  the  oscillatory  AICs  required  the  assumption  of  a 
reduced  frequency,  the  velocity 

V  =  U)b^/k^  .  (IS)  * 

The  frequency  response  problem  requires  the  addition  of  the  harmonic  gust 
forces  of  Eq.  (3)  to  the  oscillatory  applied  loads  of  Eq.  (11b).  The  total  applied 
force  becomes 

|F^|  -."(iM)  +  pb2s[C^)|hJ  +  |fJ  ,  (16) 

Substituting  Eq.  (16)  into  Eq.  (6),  and  again  introducing  the  structural  damping, 
but  here  as  its  actual  value,  permits  the  solution  for  the  frequency  response  of 
the  structural  deformation  mode 

M  .  .VW^b,.lRl[cJ|w^/W^|  (17, 

where 

[R]  =  ([I]  -  {J/{1  +  ig^))[a]  ([M]  +  pbJs[Cj^]))  .  (18) 

If  the  response  calculation  is  carried  out  for  constant  velocity,  then  the  range  of 
reduced  frequency  used  to  determine  the  oscillatory  AICs  for  motion  and  for  the 


gust  follows  from  the  range  of  frequency  £  w  <  of  interest 

w.b  /V<k  <a).,b  /V  (19) 

1  r  r  N  r  '  ' 

From  the  frequency  response  of  the  structural  deformations,  the  frequency  re¬ 
sponse  of  any  stresses  may  be  found,  and  a  basis  is  provided  for  random  loading 
fatigue  calculations. 

The  above  examples  have  shown  the  application  of  AICs  in  a  collocation 
formulation  of  various  aeroelastic  problems.  But  a  collocation  solution  is  not 
always  feasible  when  a  very  large  number  of  degrees  of  freedom  must  be  treated 
and  a  series  or  modal  solution  of  the  collocation  equations  must  be  sought.  The 
modal  formulation  may  be  derived  a  number  of  ways  (e.  g.  ,  from  the  collocation 
equations  using  the  Galerkin  method,  or  from  the  mass,  stiffness,  and  AIC  ma-* 
trices  using  the  Lagrange  method).  We  do  not  wish  to  illustrate  the  modal  form¬ 
ulations  of  the  preceding  examples,  for  these  are  perhaps  better  known  than  the 
collocation  formulations.  However,  we  do  wish  to  illustrate  the  calculation  of 
the  generalized  aerodynamic  forces  from  the  AICs  for  use  in  the  modal  approach. 
Consider  a  series  approximation  to  the  deflection  mode, 


|h|  =  rh„]  |q„|  (20) 

where  N  is  a  rectangular  matrix  of  the  assumed  modes  (e.  g.  ,  each  column 
may  be  one  of  the  free -vibration  modes),  and  the  are  their  respective  gen¬ 

eralized  coordinates.  We  may  compute  the  generalized  force  in  each  mode  by 
considering  the  virtual  work  done  in  perturbing  the  mode.  The  virtual  work  done 
in  perturbing  the  m'th  mode  is  given  by 


and  the  generalized  force  is  given  by 


=  (5W^/6q^ 
m  m  m 


6h  /6q 

m 

1 
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We  note  from  the  assumed  seriest  Eq.  (20),  that 


with  which  we  obtain 


(23) 


(24) 


We  may  extend  this  result  to  obtain  all  of  the  generalized  forces 

|q|  =  1f|  .  (25) 

Combining  Eqs.  (20)  and  (25)  first  with  Eq.  (If  leads  to  the  static  generalized 
force  s 


|q^|  =  (qS/E)[hJ'^[ChJ[hj|qn|.  .  (26) 

and  then  with  Eq.  (2)  leads  to  the  oscillatory  generalized  forces  for  motion 

1q|=  |q„|  ;  (27) 

combining  Eqs.  (3)  and  (25)  leads  to  the  oscillatory  generalized  gust  forces 

l“g|  • 

At  this  point,  it  is  obvious  that  the  AICs  offer  a  large  computational  advantage 
in  aeroelastic  analyses  by  modal  methods,  since  the  AICs  need  be  derived  only 
once,  and  the  generalized  forces  follow  from  multiplications  with  the  modal  ma¬ 
trix;  any  appraoch  that  derives  generalized  forces  directly  is  inefficient  to  the 
extent  that  the  aerodynamic  analysis  must  be  repeated  for  every  revision  of  the 
modal  set. 

The  discussion  throughout  this  paper  will  center  on  theoretical  results.  We 
cannot  discuss  experimental  correlation  of  the  theoretical  AICs,  since  sufficient 
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unsteady  measurements  are  not  available.  Any  modification  of  theoretical 
unsteady  AICs  must  be  based  on  the  limiting  case  of  the  static  AICs  and  corre¬ 
sponding  steady  wind  tunnel  measurements.  A  derivation  of  a  weighting  matrix 
that  permits  adjustment  of  the  theoretical  AICs  to  incorporate  the  results  of  ex¬ 
perimental  measurements  is  presented  in  Section  V. 


n.  AICs  FOR  MOTION 


A.  Strip  Theory 

We  consider  first  the  derivation  of  the  AICs  from  the  simplest  of  the 
aerodynamic  theories  -  -  that  of  strip  theory  in  which  the  flow  along  any  sec¬ 
tion  of  the  wing  can  be  considered  as  two-dimensional.  The  derivation  is 
taken  from  Ref.  1,  and  is  repeated  here  to  illustrate  the  three  basic  relation¬ 
ships  that  must  be  established  for  any  aerodynamic  theory  to  obtain  the  AICs; 
(1)  the  pressure -downwash  relation;  (2)  the  force -pressure  relation;  and  (3) 
the  do’vnwash -deflection  relation.  Although  these  are  the  basic  relationships, 
they  are  found  in  an  equivalent  manner  in  the  case  of  strip  theory,  since  the 
pressure  integration  has  already  been  carried  out  for  the  tabulation  of  the 
oscillatory  coefficients.  However,  this  equivalence  is  apparent,  and,  since 
the  strip  theory  approach  is  so  fundamental,  we  shall  find  it  a  satisfactory 
introductory  illustration. 


In  the  case  of  strip  theory,  the  matrix  of  AICs  appears  in  a  partitioned 
form.  For  example,  the  AICs  for  the  two-strip  wing  appear  as 


(29) 


where  the  Cj^^  are  the  AICs  for  each  strip.  Hence,  it  is  only  necessary  to 
derive  the  AICs  for  one  strip.  The  given  and  the  equivalent  force  systems  and 
geometry  are  shown  in  Fig.  1.  We  have  arbitrarily  chosen  the  forward  control 
point  at  the  quarter -chord,  the  intermediate  control  point  at  the  control-sur¬ 
face  hinge  line,  and  the  aft  control  point  at  the  trailing  edge.  The  chord  is  con¬ 
sidered  to  be  rigid.  The  oscillatory  aerodynamic  coefficients  referred  to  the 
quarter -chord  are  defined  by 
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■  Fig.  1.  The  Original  (a)  and  the  Equivalent  ( b)  Force 
Systems  and  Geometry  for  Strip  Theory 


(N.B.,  the  factor  tt  cos  A  is  used  with  subsonic  coefficients;  usually  the  factor 
4  is  used  in  its  place  with  supersonic  coefficients.)  Eq.  (30)  is  the  strip 
theory  equivalent  of  the  basic  relationship  between  pressure  and  downwash. 

The  equivalent  of  the  force- pressure  relationship  arises  from  the  transformation 
of  the  given  loads  to  the  control  point  forces.  This  load  equivalence  is  given  by 
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relationship,  is  |;iven  by 


Placing  F  in  te  rms  of 


h 

■  1 

0 

0 

ba 

^  = 

-b/d 

b/d 

0 

bp 

b  4 

1  I 

b/d 

-(b/d  +  b/c  ) 

cL 

b/ c 

a  . 

(32) 


result  with  Eq.  (2)  yields  the  strip  AIC  matrix 
(Ch)  =  TTcosA(b/ (Ay/s) 


h  by  combining  Eqs.  (30)  to  (32)  and  identifying  the 


In  the  absence  of  the  control  surface,  we  have 


2 

1  -b/d 

L  L 
h  a 

1  0 

=  ■iTcosA(b/l;>^  )  (Ay/s) 

0  b/d 

h  a 

-b/d  b/d 

(34) 


Once  the  AIC  partitions  are  obtained  for  each  strip,  the  surface  AICs  are 
formed  by  assembling  the  partitions  according  to  Eq.  (29). 

A  survey  of  the  status  of  two-dimensional  oscillatory  aerodynamic  theory 

must  begin  with  the  incompressible  solutions  of  Theodor  sen  and  of  Theodor- 

1415  16 

sen  and  Garrick  *  and  the  tabulations  of  Smilg  and  Wasserman  .  These 

solutions  are  for  the  case  of  a  rigid  airfoil  and  control  surface;  the  incom¬ 
pressible  solution  for  the  case  of  a  flexible  chord  undergoing  parabolic  changes 
in  camber  has  been  obtained  by  Spielberg  The  subsonic  solution  fer  the 
rigid  airfoil  and  control  surface  has  been  obtained  and  tabulated  by  Timman, 
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van  de  Vooren,  and  Griedanus.  ''  ,^onic  solution  has  been  given  by  Rott, 
and  the  extension  to  include  ?  v  surface  and  a  tabulation  have  been  given 

ay  Nelson,  and  Berr.nan.^°  ^^^upersonic  solution  for  a  flat  plate  airfoil  in¬ 
cluding  a  control  su*-^  been  obtained  by  Garrick  and  Rubinow^^  and 

tabulated  iiy  Hu'-'  ^  Curling.  The  zero  thick;iess  supersonic  solution 
tends  to  prr  ,  V^servative  flutter  predictions  beciuse  of  its  aft  location 
of  ae-  .  ■  ’  center.  A  supersonic  frequency  e>c.,3nsion  solution  for  the 
pressure  including  thickness  effects  but  without  a  control  surface  has 
.  given  by  VanDyke.^^  and  the  pressure  coefficient  has  been  integrated  by 
1)  ,r  ,den  and  Revell  to  obtain  the  oscillatory  coeffi  -  ents;  the  results  have  not 
b^on  tabulated.  The  solution  with  thickness  effects  approaches  second-order 
/liston  theory  at  high  Mach  number.  Third -order  p  ton  theory  including  a 
Control  surface  has  been  discussed  by  Ashley  and  Z  .rtarian;^^  this  solution 
•vill  be  reviewed  in  a  three-dimensional  context  in  s  later  section. 

•3.  Lifting  Surface  Theory 

1.  Subsonic  Lifting  Surface  Theory 


Introduction 


JiiJhe  f olio wAi 


the  aerodynamic  influence  coef^SHHtow^?^*Ti?ting  surface  Repre- 

renting  the  control  point  force  to  deflection  relation  I o  be  as  pr  evious-l v 
in  Section  I ,  ■* 


sections  we  dedeci  tl 


flit  in g 


equations  for 


-  p(ju 


hb  fej 


then  wc  shall  show  that  a  convenient  representation  s 

[cil=  K(l/kJ)  |bJ  |w  J  (3 

where  ^  pressure  integration  matrix,  J  matrix  relating 

•pressure  loading  coefficiemts  to  downwash,  and  J  relates  downwash  to 

deflections  of  the  collocation  points.  We  first  cons.vler  these  component  ma - 
rrict's  in  detail.  • 


'I'We  refer  to  methods  leading  to  an  integral  ecjuation  relating  dov/nwash 
and  the  pressure  loading  on  the  wing  represented  b\  a  series  of  functions  with 
unknown  coefficients  to  be  determined  as  the  main  f-roblem. 
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b. 


Review  of  Earlier  Solutions 


The  general  theoretical  foundation  for  the  solution  of  the  sub- 

26 

sonic  unsteady  problem  was  laid  by  work  of  Kuessner,  and,  subsequently, 

27 

by  Watkins,  Runyan>  and  Woolston.  With  the  advent  of  large  high-speed 
digital  computers,  the  development  of  numerical  methods  followed.  Among 
these  was  the  method  of  Runyan  and  Woolston,  which  was  subsequently 
developed  by  Rodden,  Hodson,  and  Revell  in  the  form  of  AICs  according  to 
the  present  formulation.  This  development  will  be  discussed  in  the  context  of 
the  general  method  to  illustrate  the  technique.  The  more  recent  work  of  Wat¬ 
kins,  Woolston,  and  Cunningham,  having  more  general  applicability,  forms 

the  basis  for  recent  developments  of  Revell,  Carrington,  Milligan,  and 
31 

Vivian.  A  development  basically  parallel  to  the  NASA  work  has  been  made 
32  33 

by  Hsu  *  and  the  M.  I.  T.  Aeroelastic  Research  Laboratory  (see  Ref.  34  for 
a  general  review  of  this  work).  All  of  these  subsonic  methods  are  based  on  the 
acceleration  potential,  or  pressure  loading  approach,  wherein  the  unknown 
pressure  distribution  on  the  planform  is  represented  by  a  series  of  functions 
whose  unknown  coefficients  are*related  to  the  known  downwash  distribution 
by  an  integral  equation,  which  for  the  harmonic  case  is  of  the  form 


w 


(x,y,0,w)  JJ  K  (x  -  g,  y  -  ti,  M,  u)  Ap  (^,  t))  d|  dT| 


(37) 


and 


M  N 

(5.  Tl)  =  ^  £  ^Pnm  ^ 


m=0  n=0 


(38) 


nm 


We  consider  the  development  oLAICs  as  illustrated  by  the  earlier  development 
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of  Rodden,  Hodson,  and  Revell  in  the  context  of  the  current  work  of  Ref.  31. 

The  derivation  of  the  AICs  from  an  arbitrary  lifting  surface  theory  is 
always  possible,  however  complex  may  be  the  details.  We  present  a  systematic 
approacl^ which  is  valid  for  cither  subsonic  or  supersonic  flow  when  the’ acceler¬ 
ation  potential  method  is  used.  This  form  is  preferred  over  the  velocity  po¬ 
tential  method  in  supersonic  cases  involving  subsonic  edges  as  will  be  discussed 
in  Section  II.  B.  3).  We  define  the  downwashes  induced  at 


a  point  set  ji|  in  the  field  to  be  linearly  related  to  the  load  distribution  ex¬ 
pressed  as  a  series  of  suitable  functions  cind  undetermined  coefficients  la  I 

J  nm  I 

(n  =  1(1)N,  m  =  l(l)M).  In  matrix  form,  we  write 


Iw  /vl  =  [d^  I  la  I 

(  1  I  L  nm]  I  nm  J 


*  (39) 


where  in  the  subsonic  case,  following  Watkins,  Woolston,  and  Cunningham 


.  b  f  r 

D*  =  °  /  dn 

nm  4tt  / 


Kq  (X.  -e,  y.  -n,  M,  k)  Ap^  (e.t,)  di 


f  (Yi  -  11)' 


,(40) 


M  r-  -5-  N 

^Pnm^P^  =  4TT{f/bQ)  (b^/b)  ^  ^  ^  JQ)  •  (41) 

m  =  0  n=0 


The  planform  geometry  is  shown  in  Fig.  Z,'  and 


i  ~  ^i^(n)  -  (b(r|)/b^)  cos  0  .  (4Z) 

Before  discussing  the  development  of  AICs,  we  pause  to  consider  the  problem 
posed  by  Eq.  (35)  which  is  the  heart  of  the  aerodynamic  problem. 


Iz 


V 

Fig.  Z.  Subsonic  Planform  Geometry 
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c.  Calculation  of  the  Induced  Downwash 

The  matrix  of  Eq.  (39)  is  handled  several  ways.  The 

Runyan  and  Woolston  method  follows  the  approach  of  replacing  the  distributed 


loads  by  a  discrete  set  of  concentrated  loads  <  AP  (x  ,7]  )>  and  related  them  to 

(  V  V  u  ) 

i  I  by  a  set  of  "downwash  factors" 
for  a  span  wise  strip  wide,  so  that 

(w./v|  =  (l/4iTpV^)[F^j^/£j^]  jAPJx^,r,^)|  .  (43) 

Now  the  concentrated  loads  are  related  by  an  integration  matrix 
the  pressure  loading  coefficients  |  |  so  that 

|a„^|  (44) 


where  from  Eq.  (41)  we  have 


Then  Eqs.  (43),  (44),  and  (39)  imply  that 


In  the  method  of  Runyan  and  Woolston,  and,  hence,  the  above  developments,  the 
downwash  factors  are  expanded  in  a  fifth-degree  polynomial  in  the  reduced  fre- 
quency  k  =  cob^/V.  Furthermore,  the  fifth-degree  term  behaves  as  k  /|3  as 
M-*  1;  whence,  the  solution  is  only  valid  when  k.Vl/p  <<1,  although  it  should 
serve  as  adequate  for  many  lowfM*  mode  instability  problems  for  Mach  numbers 
up  to,  say  0.80.  Application  of  the  ratio  test  shows  convergence  whenever 

I  kx„  |<  lOp^/3,  which,  in  general,  is  nonuniform  on  the  planform  (through  the 

^  II  2  - 

Xq  variation).  However,  for  delta  wings,  Ix^KZ,  and  k  <  5p  /3  is  a  sufficient 

condition  for  uniform,  absolute  convergence. 

The  method  of  handling  1^^^]  in  Ref.  Z9  differs  somewhat  from  that  of 
Runyan  and  Woolston,  but  leads  to  similar  quantitative  results.  An  improved 
computational  procedure  for  this  method  has  been  developed  by  Rodden,  Farkas, 
and  Oyama. 
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The  more  recent  works  of  Watkins,  Woolston,  and  Cunningham"^  and 
Hsu^^*  proceed  to  evaluate  [pj^  J  directly  by  Gaussian  quadrature  for  the 
nonsingular  regions,  and  uses  the  finite  part  concept  of  Hadamard  to  treat  the 
singular  region  (seeNRef,  36,  pp.  Z77-Z80).  The  development  of  Revell, 
Carrington,  Milligan,  and  Vivian^^  considers  piecewise  continuous  planform 
edges,  and  the  integration  must  be  interpreted  in  the  Lebesque  sense.  The 
integration  scheme  of  Ref.  30  is  essentially  valid  over  the  entire  frequency 
range,  and  the  more  accurate  integration  provides  a  more  reliable  represen¬ 
tation  of  the  effect  of  the  distributed  loading  than  does  the  replacement  loads 
method.  This  matrix  is  the  point  of  departure  for  the  remainder  of 

the  problem,  and  the  question  arises  as  to  how  many  loading  functions  and 
downwash  match  points  to  demand  for  an  adequate,  efficient  engineering  so¬ 
lution  of  the  problem. 


d.  Numerical  Solution  of  the  Downwash  Integral  Equation 

An  exact  solution  of  the  linear  boundary  value  problem  of 
lifting  surface  theory  demands  that  the  downwash  shall  match  some  prescribed 
values  everywhere  on  the  planform.  The  boundary  condition,  therefore,  requires 
N^“*-oo.  For  computational  purposes  we  settle  for  the  finite  number  of  points 
^i^^i'  ^i^  ‘  Historically,  the  early  computational  procedures  developed  by  Wat¬ 
kins,  Woolston,  and  Cunningham  have  chosen  N.  =  9,  and  16.  The  method  of 

32  33  .  .  ^ 

Hsu  '  involves  a  25-point  collocation  scheme  whose  locations  are  chosen 

to  preserve  the  total  chordwise  load.  In  order  to  solve  Eq.  (39),  it  is  clear 
that  we  need  one  equation  for  each  unknown  a  .  Thus  the  simplest  view  of 
the  problem  is  to  set  N.  =  N  x  M  and  make  as  large  as  possible  for  a  given 
digital  computing  capability.  However,  because  of  the  polynomial  form  of 
the  matrix  is  not  "well  conditioned”  for  a  finite  digit  arithmetic  system. 

A  common  remedy  to  such  problems  is  the  least  square  solution  which  over- 
determines  the  problem  (see  Ref.  37,  pp.  149-164).  Thus,  we  shall  choose  a 
relatively  small  value  of  N  x  M  which  is  sufficient  to  describe  our  load  distri¬ 
bution,  but  endeavor  to  increase  N.  as  far  as  possible.  We  thus  find  a  non- 

finalitv  of  the  (a  !,  but  for  each  combination  N  &  M  we  find  a  convergence  of 
/  I  nm)  -  - - 


the  {a  }  as  N.^oo. 
I  nm)  1 


15 


1% 
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© 


We  thus  expect  to  find  an  optimum  computational  solution  by  choosing 
(Nx  M)<N^  to  reduce  the  size  of  the  associated  algebraic  problem.  This  point 
is  clear  when  we  consider  the  general  solution  of  Eq.  (39)  for  the  case  N^>(N 
X  M).  Here  we  havq 


a  =  A' 

nm  I  nrn 


w./V 

1 


(47a) 


where  from  Eq.  (39) 


nm 


nm 


-1 


(47b) 


where  JA  Sjj  is  the  least  squares  diagonal  weighting  matrix  whose  elements 
consist  of  the  incremental  areas  surrounding  the  control  points.  We  note  that 


D^mJ  complex  for  the  oscillatory  aerodynamic  problem,  and  it  can  be  con¬ 
sidered  as  the  Fourier  transform  of  the  general  time  dependent  downwash. 
Noting  the  ranks  of  the  matrices  carefully,  we  know  that  the  solution  of  Eq. 

(47)  exists  only  if  N^XN  x  M).  Physically  speaking  we  are  approximating  the 
point  boundary  condition  in  a  least  square  error  sense  for  any  given  N  x  M 
combination  of  loading  functions.  We  shall  now  discuss  the  loading  functions 
briefly. 


e.  Choice  of  Loading  Functions 

(1)  Chordwise  Functions:  'I’he  most  widely  used  loading 
functions  satisfy  the  Kutta  condition  of  pressure  continuity  at  the  trailing  edge, 
and  have  the  classical,  tliin  airfoil  square  root  singularity  at  tlie  leading  edge, 
lu' terms  of  the  trigonometric  variable,  P  ,  defined  by  Eq.  (43)  the 
Eq.  (41)  arc  of  the  form 

(P)  -  rln  (P/Z)  (48a) 


^  (B)  -  (4/Z^”)  sin  (nP), 


11  ■>  1 


(48b) 


hi  the  case  of  a  trailing  edge  flap,  we  have  a  discontinuity  in  downwash, 
and  Jtfsu'^'^  has  shown  that  a  logarithmically  singular  loading  term  arises  of  the 
form 


A  w„ 

I  (k) 
c 


M 


Z 


I  -  cx  I  \/l-x^ 

1  -  cx  -  yjx  - 


(t  -  ^ 


{4'>) 


If) 


Z  (  f  /ll^) 


for  a  control  surface  covering  c<  jc  <  1,  and  the  wing  chord  is  defined  by 

-1<  x<  1. 

(2)  Spanwise  Loadings;  Interference  Effects;  For  the 

treatment  of  interference  effects  it  is  convenient  to  modify  the  spanwise  load¬ 
ing  function  given  by  Eq.  (45),  for  clearly  at  the  wing  root,  no  combination  of 
functions  of  the  type  will  alter  the  loading  for  Falkner^  has 

contrived  the  so-called  P  functions  for  this  purpose.  However,  for  the  pur¬ 
poses  of  thisdiscussion,  any  convenient  analytic  function  which  is  even  in  T)-q^, 
where  is  the  center  of  a  spanwise  interference  effect,  w^uld  be  sufficient. 

We  might  choose,  for  example 

8m  \/l  -7}^  cos  mn  (50) 

where  M>m,  or 


Having  described  the  classes  of  convenient  subsonic  loading  functions  which 
satisfy  the  boundary  condition,  we  consider  the  problem  of  downwash  specifi¬ 
cation. 


f.  Specification  of  the  Downwash;  Interference  Considerations 


The  planform  boundary  condition  for  the  general  aerodynam¬ 
ic  problem  is  to  match  the  normal  component  of  velocity  of  the  fluid  relative  to 
the  surface 


(52) 


and  for  the  linearized  lifting  surface  problem,  this  amounts  to  evaluating  w- 
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on  the  plane  z  =  0,  (cf.  ,  for  example.  Miles  ).  Hence,  Eq.  (39)  expresses  the 
fluid  downwash  velocities,  w  induced  l)y  the  wing  loading  at  any  point,  (x^, 

V.)  in  the  plane  z  =  0.  The  body  vertical  velocity,  at  (x^,  y^),  therefore, 
must  equal  the  sum  of  all  the  induced  fluid  downwash  velocities  from  each  of 
disturi^ancc  fields  (of  which,  the  wing-induced  downwash  is  the  major 


component).  Then 


The  body  downwash'for  an  arbitrary  surface  deformation  pattern  Z  =  h(x,  y,  t) 
is  to  the  first  order  of  approximation, 


w  =  (V  —  + 
bi  \  Sx  a 


X,  y=x..y. 


For  the  purpose  of  the  collocation  method  as  discussed  in  the  previous  sections, 
we  must  relate  deflections  at  a  discrete  grid  point  set  j  h^  j 

(  j  =  1  (l)Nj).  To  this  end  we  define  J  such  that 


[w^./Vj  =  (l/bo)  [W..)  Ih.J  _ 

The  problem  of  computing  j  requires  some  type  of  analytic  inter¬ 

polation  of  the  deflections.  For  example,  let 

h  (x,y)  =  ^  Y  ^rs  “rs  . 

r  s 

At  our  control  point  grid  we  have 

IM  =  K.|1--I  ■ 


h .  =  IF 

J  ’ 


Using  a  local  chordwdse  interpolation  only,  we  obtain  for  y  =  y^  ,  x  = 

|'\i  <yi>  I  =  ["i  <^j>l  l^’r  <yi)|  .  J  =  I  (I)  N.  (yj) 

Then  by  least  squares 

K''i'i= Kiy  i»ii  hJiiM 


■ii  r'j 


wiicre  the  least  squares  diagonal  weighting  matrix  is  composed  of  the  incremental 
chord  lengthsZ^Xj,  'Fhen  EcjS.  (59)  and  (55)  imply 


lW;.(yi: 


+  ikH 


A  general  discussion  of  this  problem  would  lead  us  too  far  afield,  but  it 
is  an  important  problem.  We  merely  remark  here  that  a  local  chord  or  interpo¬ 
lation  "in  the  small"  is  most  advisable  with  considerable  "overdetermination" 
of  the  j  j  being  a  necessary  requisite  for  reliable  results.  Interpolation  of 
the  entire  surface  ("in  the  large")  can  often  be  done  but  frequently  leads  to 
nearly  singular  matrices. 

To  consider  an  interference  problem,  we  specialize  Eq.  (53)  to  the  case 
of  a  wing  and  a  fuselage  and  write  using  Eq.  (39), 


(61a) 


„i\V  1  I  W 

D  a 

nm  j  I  nm 


j  I 


(61b) 


Then  our  solution  for  Eq.  (47)  requires  that  we  define 


W/..  I  W...  WF... 

iw.  /V  =iw,./V  -  w.  /V 

II  i  bi  I  i 


(62a) 


W  \V 

(1/1, q)  ]  1,' 


/  V 


(62b) 


wjicrcjw.  /v|  might  be  evaluated  from  some  separate  analysis  and  related  to 

the  local  fuselage  downwash  and  deflections.  For  example,  by  slender  body 
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theory  the  fuselage  field  perturbation  potential  is  (cf.  ,  Miles  ),  noting  sin  P  = 
z/r  for  a  circular  fusedage 

(x,t;  r,P)  =  =-  (l/r)w|^  (x,t)  Rp(x)  sin  0  ,  (6 

whence,  the  downwash  in  the  field  is 


WF  F 

w.  (x^;  r.,  0.,  t)  =  3  /a  /. 


(64a) 
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(64b) 


2  7 

For  z  «  r 
1  1 


=  (l/rf)  wj^  (x..t)  (x.)  (1  -  2zf/rf) 


WF  F  ,  02  /  X/  2 


Further  if  the  fuselage  downwash  is  related  to  the  deflections  by 


(x„t)  -  (1/b,)  W"  ht 


1  •  r  2  2 

we  have  again  for  z  .  «  r.  , 


|wf''^-(x.,t;y.,0)/v|  =  (1/b^)  lR^(x.)/rfl  \wf.\  |h[| 

Then  substituting  Eq.  (66)  into  Eq.  (62b)  gives  an  aerodynanaic  interference 
coupling  between  the  fuselage  and  wing  motion 

|wf/v|  =  (l/b^)(lw;;!  |h]"|  -  [R^(x.)/rf]  |wj|  |h[|)  ( 


which  may  be  written  in  the  partitioned  form 


rr/vi  -  K*  I  "rl  -> 


(67b) 


r  w 1  f  WF 1 

The  matrices  J  jj  J  obtained  by  means  of  deflection  interpo- 

lationas  discussed  previously.  We  next  consider  development  of  the  matrix 
relating  the  forces  j ^ j  j  loading  coefficients,  | 


Relation  Between  Control  Point  Forces  and  Pressure  Loading 


Coefficients 


This  is  perhaps  the  most  intuitively  simple  aspect  of  the 
method,  and  in  the  limit  of  an  infinite  number  of  force  points  gives  |cf.  ,  Eq.  (41)j 

Fj  ~  APj  ASj  ~ /)V^  /  ^  (e.)  g^(Ti.)  As^  .  (6h) 

However,  what  we  are  interested  in  is  a  finite  lumping  of  the  aerodynamic 
loads  which,  in  some  sense,  preserves  the  most  important  properties  of  the 
distributed  load  and  must  be  placed  at  points  chosen  arbitrarily  by  the  dynam¬ 
ics  engineer  (for  reasons  of  inertia  and  stiffness  data  availability).  We 

formalize  the  definition  of  B  by  the  relation 

I  nm  J  ^ 
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p.l  = 


irp 


l^L] 


nm 


(69) 


The  simplest  treatment,  and  one  which  is  especially  appropriate  for  a 
large  number  of  chordwise  collocation  stations,  is  to  write 

hL|  =  SS  dT,.  j  -  1(1)N.  .  (70) 

As. 

J 

In  the  method  of  Ref.  31  for  the  cases  where  41N.(y)  <  5,  the  strip  lift, 
pitching  moment,  and  control  surface  hinge  moment  are  matched,  and,  rather 
than  matching  higher  order  moments,  an  alternative  method  called  the  "mini¬ 
mum  vector  length  constraint"  is  used.  Although  the  method  of  matching 
higher  order  moments  is  necessary  to  yield  the  correct  generalized  forces 
for  an  arbitrary  chordwise  deformation  mode,  it  does  not  lead  to  a  satisfactory 
estimation  of  the  chordwise  load  distribution  for  the  purpose  of  stress  analysis 
when  only  a  small  |e.  g.  ,  Nj(y)  “  ^  of  chordwise  control  points  are 

chosen.  The  minimum  vector  length  constraint  adds  to  the  basic  conditions 
of  lift,  pitching  moment,  and  hinge  moment,  the  requirement  that  the  sum  of 
the  squares  of  the  control  point  forces  be  a  minimum.  We  estimate  that  the 
procedure  defined  by  Eq.  (70)  is  satisfactory  for  ^^(y)  of  the  order  of  six  or 
more  per  half  wave  of  chordwise  deformation,  and  the  minimum  vector  length 
constraint  approach  has  its  utility  in  nearly  rigid  chordwise  deformation  prob¬ 
lems.  We  add  that  one  must  keep  a  balanced  view  in  regard  to  aerodynamic 
lumping  versus  the  inertial  lumping  for  the  vibration  analysis,  and  with  due 
consideration  for  the  quality  of  the  structural  influence  coefficients  for  the 
aeroelastic  analysis. 


h.  The  Aerodynamic  Influence  Coefficients 

We  are  now  in  a  position  to  assemble  the  component  matrices 
to  oli^ain  the  desirer!  AICs.  The  fundamental  relationships  as  described  in  Sec¬ 
tion  II  A  are  formed  in:  Eq.  (47a),  the  pressure -downwash  relation,  Eq.  (69), 
tlie  force -pressure  relation,  and  Eq.  (55),  the  downwa sh -deflection  relation. 
Combining  these  so  as  to  relate  the  forces  and  deflections  leads  to 
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(71) 


|Fj|  |w,.|  |h.|  . 

Identifying  Eqs.  (35)  and  (72)  leads  to  the  AICs 

[cj'  =7r(i/k5)(Wbo)  [a|^^|  [W..j  .  (72) 

We  have  elaborated  on  the  formulation  for  the  subsonic  case  because  it 
contains  most  of  the  possible  complexities  of  applying  lifting  surface  theory, 
and  we  have  demonstrated  the  feasibility  of  the  collocation  format  of  the  AICs. 

We  next  consider  lifting  surface  theory  in  other  Mach  number  regimes. 

2.  Transonic  Lifting  Surface  Theory 

a.  Introduction 

This  regime  is  the  most  important  of  all  the  Mach  number 
ranges  in  aeroela stic ity,  even  for  vehicles  whose  design  Mach  numbers  are  in 
the  supersonic  range.  In  fact  even  for  a  hypersonic  vehicle,  there  will  exist 
two  critical  regimes  controlling  stiffness: 

(1)  A  transonic  design  point,  due  to  the  fact  that  aerodynam¬ 
ic  derivatives  reach  their  maximum  in  this  regime. 

(2)  A  hypersonic  design  point,  due  to  the  forward  shift  of 
aerodynamic  center  due  to  thickness  effects,  and  due  to  temperature  effects  on 
stiffnes  s. 

Perhaps  the  greatest  source  of  difficulty  in  the  transonic  regime  on  a 

hypersonic  vehicle  is  the  type  of  configuration  (e.  g.  ,  blunt  leading  edges  and 

3  8 

very  low  aspect  ratios,  such  that  the  planar  surface  concept  breaks  down, 
and  we  must  proceed  with  a  "wing-body"  approach  from  the  outset).  However, 
for  all  vehicles  for  which  aerodynamic  heating  is  not  crucial  in  dictating  wing 
design  (i.  e.  ,  approximately  up  to  M  =  3.  0),  we  shall  face  the  problem  of  low 
aspect  ratio,  very  thin  wings,  since  wave  drag  is  an  important  consideration 
in  design.  For  such  vehicles  the  transonic  problem  is  essentially  three  di- 

39 

mensional  and  lifting  surface  theory  is  the  only  reasonable  basis  for 
aeroelastic  loads  analysis.  • 
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b. 


Linearization;  Shock  Wave  Effects 


It  is  generally  appreciated  that  shock  waves  and  nonlinearity 

are  the  essential  features  of  the  steady  transonic  flow  past  a  body  or  surface  of 
40'  38  39 

finite  thickness,  but  it  has  been  shown  *  that  there  exists  a  valid  range 

4 1 

\  for  linearization  for  the  unsteady  regime,  for  low  aspect  ratio  wings  defined  by 
\  (see  Fig.  3) 

k»  I  at  log  (at  (73) 

where  O  =  l/Zh^,  €  =  thickness  ratio,  and  k  =  ^b^/V  is  the  reduced  frequency. 

(For  example,  for  a  60-degree  swept  delta  wing  with  f  =  0.  03,  O  =  1/^37  we 

require  k  »  0.  0297).  In  this  event,  the  appropriate  linearized  equation  for  the 
,  ,  .  38,39 

unsteady  motion  is 


*^1  +  -  2M^a!  -  «S  =0  .  (74) 

yy  zz  xt 


Fig.  3.  Planform  Geometry  for  Transonic  Regime 


The  presence  of  shock  waves  has  an  important  effect  in  modifying  the 

42, 

boundary  conditions.  Eckhaus  has  recently  developed  a  solution  of  the  ideal¬ 
ized  problem  of  an  unsteady  two-dimensional  flow  with  a  uniform  supersonic 
zone  separated  from  a  uniform  subsonic  zone  by  a  normal  shock  wave  perpen¬ 
dicular  to  the  airfoil  surface  (see  Fig.  4a).  This  actually  is  a  fairly  realistic 
model  of  the  subsonic-transonic  regime  (see  Fig.  4b). 


1  * 

su 

M|  >  1  1 

M2  <1 

—  V. 

1 

1 

-  - 

j 

1  XSL 

SHOCK  TRACE  IN 
PLANE  y=  CONST 


SONIC  SURFACE  TRACE 
IN  PLANE  y=  CONST 


(a)  IDEALIZED  FLOW 


b)  SUBSONIC-TRANSONIC  FLOW 


Fig.  4.  Transonic  Flow  with  Shock  Waves 


Eckhaus'  model  has  already  proven  to  be  helpful  in  explaining  many  of  the 
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features  of  control  surface  buzz  *  wherein  a  disturbance  due  to  control  sur¬ 
face  oscillation  sends  upstream  a  wave  which  interacts  with  the  shock  in  a  man¬ 
ner  such  that,  as  always,  the  Rankine -Hugoniot  relations  must  be  satisfied. 

For  the  unsteady  problem  these  provide  a  set  of  boundary  conditions  relating  the 

39  42  43 

unsteady  perturbation  potential  derivatives  on  either  side  of  the  shock.  '  ' 

Eckhaus  shows  that  the  most  important  first-order  effect  of  the  shock  is  to  make 
the  airfoil  appear  "semi-infinite"  (with  respect  to  an  origin  fixed  in  the  steady 
state  shock  wave  position).  Part  of  the  solution  yields  the  amplitude  of  the 
shock  wave  oscillation  iti  relation  to  the  surface  motion  which  originated  the 
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disturbance.  This  is  also  of  interest  in  connection  with  problems  of  shock  wave 

interaction  with  turbulence  which  are  of  importance  in  the  growing  field  of  aero- 

44 

dynamic  noise  theory  pioneered  by  LighthilL  For  a  review  of  these  turbu- 
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lence  interaction  problems  the  reader  is  referred  to  Chu  and  Kovasznay. 

The  Eckhaus  method  could  provide  a  useful  basis  for  the  improvement  of 
lifting  surface  theory  as  a  prelude  to  development  of  AICs  if: 

(1)  The  steady  state  flow  field  about  the  wing  body  is  given 
or  can  be  estimated  (say,  for  a  symmetrical  body); 

(2)  The  shape  and  location  of  the  shock  surface  is  given  by 
an  equation  of  the  form  F^(x,  y,  z)  =  0; 

(3)  The  Eckhaus  solution  is  extended  to  treat  the  case  of 
upstream  disturbances  and  their  transmission  through  the  shock  wave; 

(4)  The  solution  is  formulated  in  terms  of  acceleration  po¬ 
tential  to  avoid  integration  over  the  wake,  in  order  to  be  more  compatible  with 
a  unified  kernel  function  approach. 

c.  Vorticity  Effects 

Wc  note  finally  that  linear  methods  do  not  preclude  considera- 

43 

tion  of  vorticity  effects,  for  we  may  write,  following  Sears, 


lere^x  v^  0,  but  where 


=  i  V  +  +  Vj 


Vj|/V  «  1  ;  I  Vq^/V 


leading  to  equations  of  the  form  (specialized  to  steady  flow  here  for  illustration 
of  the  influence  of  steady  flow  nonuniformity) 


[  1  -  (y,  7.)  j  - 


2M  (y,  7,) 


2M  (y,  z) 


M(y,  7.)  M  (y,  z) 


P'  +  P:.„  +  P’  =  0  (77) 


46 

where  p*  is  the  disturbance  pressuri*.  Lighlhill  has  also  introduced  a  modi¬ 
fied  potential  function  V  such  that  (u*,  v*,  w')is  the  disturbance  velocity  field) 


^  -  P'/PqV 
Ox  ^ 


(78a) 


(78b) 


and  sliows  that  for  the  two-dimensional  case, 


(1  -  M  )  y  +  y  =0  (79) 
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This  type  of  consideration  has  been  used  by  Miles  in  studying  boundary  layer 
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effects  on  panel  flutter  and  by  Lighthill  in  shock  diffraction  effects  (e.  g.  ,  an 
atomic  blast  passing  a  corner  of  an  airfoil). 

The  above  remarks  on  vorticity  have  been  included  here  as  appropriate  to 
the  context  of  departures  from  linear  potential  theory,  but  actually,  the  vorticity 
problem  per  se  can  be  safely  ignored  for  transonic  aeroelastic  problems  (ex¬ 
cept  for  the  proper  handling  of  the  shock  boundary  conditions).  Thus  on  either  ^ 
side  of  the  shock  the  velocity  may  be  described  by 


^  ^1  y»  t)  , 


X  <  X 


(80a) 


-  i  -1  ^ 

and  the  isentropic  pressure  equation  may  be  used, 

Ap,  =  p,  Vj  +  p,  4^ 


Ap2  =  p,  4^  + 


Pz  '^z 


Linearized  Transonic  Lifting  Surface  Theory' 


(80b) 


(81b) 


If  we  consider  only  the  first-order  effects  of  the  shock  bound- 
.iry  conditions,  the  downstream  disturbances  are  excluded  and  the  linearized 

=!-The  entropy  is  p  roportiojia  1  to  j  (  y- 1 )  lyi  1  -6  ^(V  - 1  )J  j^p/^y  and  the 

ratio  of  reflected  to  incoming  wave  pressure  amplitudes  is  D^/D^  =  (1-M)  |3- 

(I  +  M)j  /(l+M)‘^[3-(I  -M)|  (Sr<-  Hof.  AZ,  pp.  Zi-ZS). 
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transonic  problem  is  obtained.  This  has  been  thoroughly  treated  by  Landahl. 

He  shows  that  a  similarity  property  exists  for  solutions  of  Eq.  (74)  described  by 
the  relation 

4^  (x,  y»  z,  t;  CT,  k,  M)  =  (1/M)  4^  (x,  M^,  M^,  t;  Ma,  k,  1)  .  (82) 
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Then  the  solution  need  be  carried  out  only  for  M  =  1.  Ashley  has  recently 
proposed  a  "transonic  box  method,  "  using  Landahl's  result 


Ashley  suggests  that  a  fixed,  square  box  grid  be  employed,  extending  out  some 
finite  distance  past  the  wing  tip,  and  that  the  diaphragm  downwashes  be  de¬ 
termined  algebraically  in  the  same  manner  as  in  the  supersonic  application  for 
the  subsonic  leading  edge  case.  (See  Section  II.  B.  3  for  a  discussion  of  this 
problem.  )  In  this  case,  the  fixed  grid  would  make  possible  an  AIC  development 
without  the  inconvenience  of  the  "floating"  Mach  box  grid  system  referred  to  in 
Section  II.  B.  3.  The  above  approach  would  constitute  a  practical  method  by 
which  Eq.  (83)  may  be  used  for  general  configurations. 

A  second  approach  to  the  linear  transonic  lifting  surface  problem  is  to 

regard  it  as  a  limiting  case  of  the  subsonic  problem.  Watkins,  Runyan,  and 
Z7 

Woolston  have  given  the  kernel  for  M  =  1,  and  calculations  have  been  car¬ 
ried  out  by  the  kernel  function  procedure.  Landahl  gives  some  comparisons 
with  his  work. 


We  finally  consider  some  limiting  cases.  For  low  aspect  ratio  and  re¬ 
duced  frequency,  we  have  the  slender-wing  theory,  for  which  the  potential  satis-  | 

fies  Laplace's  equation  in  planes  transverse  to  the  body  axis.  This  solution  has 
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been  developed  in  the  form  of  AIC  s  by  Rodden  and  Revell,  and  will  be  discussed 

at  length  in  the  following  section.  The  other  limiting  case  is  the  high-frequency 

50 

slender-wing  theory  developed  by  Merbt  and  Landahl,  in  which  the  two- 
dimensional  wave  equation  in  the  cross-flow  plane  is  solved  in  terms  of  Mathieu 
functions,  after  transformation  to  elliptic  cylinder  coordinates.  (See  Ref.  38, 

Ch.  9  for  a  review  of  this  prol)lem.  ) 
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where 
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Fig.  6.  Geometry  of  an  Elemental  Area 

The  fore  and  aft  limits  on  are  denoted  by  and  respectively,  and  the 

inboard  and  outboard  limits  on  y.  are  denoted  by  y..  and  y.  ,  respectively. 

J  J  ^  J 

The  limits  x.^  x.  ,  and  y..  are  considered  constant.  The  limit  y,  may  either 

if  la  ^ )i  ^jo  ' 

be  constant  as  in  an  interior  region,  or  be  variable  as  at  the  edge.  The  force 
assumed  to  act  at  the  point  (x.  ,  y.)  is  found  by  integrating  the  lifting  pressure 
over  the  elemental  area. 


F(x.  ,y^.  ,t)  =  J  J  Ap(x,  y,  t)  dxdy 
A  S.j 

2  f'j-*’ 

=  c  J  J  Ap(^,  n,  t)  di  dn 

where  we  denote  by  possibility  of  the  limit  varying.  The  force 

becomes 


(91b) 


F(x.  ,  y  t)  =  2p  c^ 

J 


^ia 


^  (a  f*  +  a*  f  )  I 
9^  mn  n  mn  m  n 


+  (a'  f’  +  a"  f  )  I  >  dr\di 
mn  m  mn  m  n( 
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(92) 


Since  o)  vanishes  on  the  edge,  we  may  interchange  the  order  of  differen¬ 

tiation  with  respect  to  |  and  integration  with  respect  to  r|.  Then,  if  we  define 


r  (e^,  r|)  =  (  I  (Ti,o)  d  ri 

n  J  ^ 


in  which  wc  emphasize  the  functional  dependence  cr  =  O  and 


(93) 


the  force  becomes 


V  -  V 


(94) 


F(x..y.,t)  =  2o  V*-  c; 


‘"II  fVe 


a  f  +  a'  f  )  J 
mn  m  mn  m 


+  (a'  f  +  a"  f  )  J-’  (?)  I 
mn  m  mn  m  n  j 


d!^ 


Finally,  if  we  define 


i(^)] 


(95) 


L 


mn 


and 


-[ 

(?.  r,) 


;,(?)  J„(?,9)d? 


K-*  (?)  =  K  (?.  n.  )  -  K  (?,  T|. J 

nm  mn  jo  mn  ji 


L-'  (C)  =  L  (e,  r|.  )  -  F  (f,  Ti;-) 

nm  mn  jn  mn  J1 


(96) 

(97) 

(98) 

(99) 


we  obtain  the  final  expressitm  for  the  local  force 

F(x.,y.,l)  =  ZoV^c^yy  jfa  f  (?.  )  I  a'  f  (?.  )'1  jj  (?.  ) 

-  fa  f  (?  .,)  I  a'  f  (?.,)1  J-'  (?.,.) 

1_  mn  m  if  mn  m  if _J  n  il 


+  a'  L-"  (^.  )  -  L'’ 

n'ln  mn  la  mn 


i  a"  [kJ  (f.  )  -  (I^  .j||  (100) 

il  J  mn  [_  mn  la  mn  ifJJ 
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which  can  be  abbreviated  to 


F(x.,yj.t)  =  2pV^  ^ 


a  + 


N'j 


+  K'j  a"  ) 


mn  mn  mn  mn  mn  mn 


where 


M  =  f  (P  .J  J  ri..)  +  f  (p.  )  J  (p.  ,  Ti.  ) 

nin  m  ^  if  n  "if  'ji  m  la  n  ■  la  'jo 

-  Ff  (?.J  J  Ti.  )  +  f  (?.  )  J  (?.  ,  •n..n 

1_  m  ■  if  n  if  JO  m  la  n  ^la  ji  J 


”10'  <  "ji']  * 

^mn^^ia’  ’^jo'  ‘  j^^mn^^if’  ^mn'^ia’ 


"ji>  ^mn'^ia-  '  [^iin'^if’  "jo'  +  ^mn'^a'  "ji>] 

We  see  from  the  last  three  relationships  that  the  matrix  elements  depend 
three  quantities 


M 

mn 


(f.  Tl)  =  f;^(f:)  ,1) 
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(101) 

(102) 

(103) 

(104) 

on  the 

(105) 


and 


N  (e.Ti)  =  f  (i)j  (e.T))  +  L  {i,r\) 


(106) 


K 

mn 


(e.n) 


(107) 


eval  iated  at  the  four  corners  of  each  elementary  area,  viz.  ,  the  sum  of  the 
forward  inboard  and  aft  outboard  values,  minus  the  forward  outboard  and  aft 
inboard  values. 

Before  proceeding  with  the  derivation  of  the  AICs  let  us  illustrate  the 
calculation  of  the  various  integrals  by  computing  the  lift  on  the  rigid  plunging 
and  pitching  delta  wing.  For  pitching  and  plunging  we  have 


h(x,  y,t)  =  h^lt)  +xa(t) 


{108a) 


or  m  =  0,  1,  n  =  0,  and  a^^  =  ^10  "  ^0  ~  ®0  =  ^ 

It  is  first  necessary  to  evaluate 


(108b) 


(109) 


The  first  integral  is  evaluated  from  tlie  substitution  \  -  G  cos  0  and  the  well 
known  formula  of  Glauert, 


and  for  gQ(M  =  1  leads  to 


cos  rO  do 
cos  0  -  cos  (j) 


sin  r  (|) 
sin  (j) 


(110) 


(111) 
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The  second  integration  is  routine,  and  we  obtain 


=  “Va^  -  71^ 


(112) 


The  integration  of  1  Icac^  to  d  (^,  7/). 


-  (CT^/2)  sin'‘(7,/a) 


(113) 


In  evaluating  tlic  integrals  ^qq»  it  is  necessary  to  consider 

whetlier  the  integration  is  carried  out  in  an  interior  or  an  edge  region,  i.  e.  , 
whetlior  1]  =  constant  or  7/ =  7/(^)  during  the  ^  integration.  For  a  delta  wing  of 
semiapex  angle  6  ilic  semispan  is  given  by  CT  =  ^  tan  6,  so  that  along  the  edge 
n(^)  =  Cr  =  ^  tan  6.  Carrying  out  the  integration  in  the  interior  regions,  and 
letting  //=  //  ctn  6  in  tlie  interior  regions,  yields 


Jo(irO=  const)  d^ 


(114a) 


=  -tan  d 


-  7/^  +  (4^/6)  sin"\ri/^) 


-  (17^/6)  log  (e+M^  -  rfy 


(114b) 


Kio(?»7/)  =  =  const)  d^  (115a) 

=  -tan^^  6[(l/24)(5^^  T]  -  +  (4^/8)  sin"  Vt?/^)]  (115b) 

=  |^q(?)  Jq(^»  7/  =  const)  d^  (116a) 


Ljq(^,7/)  =  j(^)  Jg(^»  7/  =  const)  d^ 


( 1  1 6a) 


(116b) 


(117a) 


=  K„nU.»/) 


(1  17l>) 
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Carrying  out  the  integration  in  the  edge  region  yields 


Ko„(e.n) 

=  J  =  &  tan  6)d| 

(118a) 

=  -  (tt/12)^^  tan^  6 

(118b) 

=  J  1  tan  6)d^ 

(119a) 

=  -  (it/  16)^^  tan^  6 

(119b) 

=J f*  (e)Jo(e»Ti  =  e  tan  6)d^ 

(120a) 

=  0 

(120b) 

LjQ(5,r|) 

=  J  =  ^tan6)d^ 

(121a) 

■  ^00 

(121b) 

Considering  the  right  side  of  the  wing  as  a  single  elemental  area  whose  four 

corners  are  the  apex  (taken  twice),  the  center  line  at  the  trailing  edge  and  the 

wing  tip,  we  may  evaluate  the  coefficients  M  ,  N  ,  and  K  at  the  corners 
^  mnmn  mn 

from  which  we  may  obtain  the  net  coefficients  and  .  For  both 

^  mn  mn  mn 

modes  considered,  all  three  quantities  vanish  along  the  center  line  and  the  net 
coefficients  become  the  values  computed  at  the  tip.  We  therefore  obtain 


Ko  = 

(DJo 

(1. 

tan  6) 

(122a) 

= 

0 

(122b) 

1 

(i)Jo 

(1. 

tan  6) 

(123a) 

(tt/4) 

tan‘ 

^6 

(123b) 
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(124a) 


^00  ^  fQ(l)  (1,  tan  6)  +  (1,  tan  6) 


=  -  (tt/4)  tan  6 


(124b) 


Nio  =  (1)  (1,  tan  6)  +  Lj^d.tan  6) 


(125a) 


=  -  (tt/ 3)  tan 


(125b) 


K  = 
^00 


=  -  (tt/  12)  tan  6 


(126) 


=  -  (tt/  16)  tan  6 


(127) 


Hence,  from  Eq.  (101)  the  total  vertical  force  becomes 


F(t)  =  2/>V^  a  +  n'-*  (c/V)a  +  (c/V)^a  )  (128a) 

/  ^  mn  mn  mn  mn  mn  rnn^ 


=  -2^V^c^  tan^  6L(iT/4)(c/V)(li^/c)  +  (n/ 12)(c/V)^  (h^/c) 


+  (TT/4)a+(TT/3)(c/V)d'+(7T/l6)(c/V)  n]  .  (128  lo) 


In  the  case  of  harmonic  pitching  about  station  x  =  ^qC,  we  let  ~ 

we  oljtain  the  lift  coefficient 


=  -(ttAR/Z)  o[l  4  ik^(4/3  -  ^^)  + 


(129) 


where  =  F/(l/4)p  V  c  tan  6  and  AR  =  4  tan  6  and  the  reduced  frequency 

is  based  on  the  root  chord,  =  coc/V.  This  result  agrees  with  that  given  by 
Miles  [Ref.  38,  Eq.(ll),  p.  357|.  The  above  digression  has  illustrated  the  evalu¬ 
ation  of  the  various  integrals.  A  complete  evaluation  of  the  integrals  for  the 
ranges  m  =  0(1)4  and  n  =  0(2)8  f(M-  the  symmetrical  case,  and  n  =  1(2)9  for 
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the  antisymmetrical  casei  is  given  in  Ref.  49»  which  provides  sufficient  data  for 
the  analysis  of  15  control  points,  e.  g.  ,  for  the  planform  divided  as  shown  in 
Fig.  7. 


To  complete  the  derivation  of  theAICs  the  matrix  counterpart  of  Eq.  (lOl)is 

|F(x.,y.,t)}  =  2pV^c^  ([M’j  ]  fa  1+  (c/V)|N*j  ] /a  \ 

V.  '  1  \[;  mn'  \  mnj  ‘  mm  \  mnj 

+  (c/V)^[K^j  ]  |a  "A  (130) 

where  the  dependence  of  on  real  time  is  now  implied.  The  matrix  counter¬ 

part  of  Eq.  (86)  is  written 


mn 


h(x.,y.,t)  =  c[H^'^  ]  /a  )  (131) 

'  1  *■  mn'  t.  mnj 

where  =  f  (£.)  g  The  solution  for  /a  \  in  terms  of  (  h}  may 

mn  m  ^n  \  mnj  V.  J 

follow  directly  from  inverting  Eq.  (120)  but  for  a  number  of  numerical  reasons 
a  least  squares  solution  is  preferable.  The  first  is  that  a  large  matrix  of  poly¬ 
nomials,  e.  g.  ,  tends  to  be  singular.  The  second  is  that  a 

Lagrangian  surface  fit,  as  we  originally  proposed,  tends  to  be  poorly  behaved 
between  collocation  points.  Finally,  it  is  often  necessary  to  consider  a  larger 
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number  of  structural  degrees  of  freedom  in  a  vibration  analysis  than  it  is 
necessary  to  consider  aerodynamic  degrees  of  freedom  in  a  flutter  analysis. 

For  these  reasons  we  consider  Eq.  (131)  to  be  over-determined,  i.  e.  ,  the  set 
|hj  is  larger  than  the  solution  to  Eq.  (131)  only  exists  in  the 

statistical  sense.  This  least  squares  solution  is 

/a  >  =  (l/c)([H’j  ]’•' [AS..]  [H'j  ]*[AS..](h)  (132) 

\  mnj  V  mn^  ij'  mnV  •-  mn^  *■  J  ' 

where  the  diagonal  weighting  matrix  of  elemental  areas,  [AS..]  ,  arises  from 

^  J 

minimizing  the  mean  square  error  over  the  surface  area.  The  oscillatory  AICs 
follow  from  substituting  Eq.  (13Z)  into  Eq.  (130),  assuming  harmonic  motion,  and 
comparing  the  result  with  Eq.  (2).  We  obtain 


=  (2/k^)  (c/s)LB'j  ]  ]*  [AS..][H*-'  ]*[AS..] 

r  mn^  mn^  ir**  mnV  '•  mm  ^  ir 


where 


] 

mn^ 


]  +  ik  (c/b  )[N*j  1  -  k^(c^/b^)[K^j  1 
mn^  r  r  mm  r'  r  ^  mn-* 


and  the  choice  of  reference  semichord  remains  arbitrary. 

The  theoretical  limitations  on  the  applicability  of  slender-wing  theory  are 
38 

outlined  by  Miles.  The  Mach  number  limitation  is  given  by 


«  (AR) 


(135) 


and  the  reduced  frequency  limitation  is  given  by 

k  (c/b  )  «  (MAR)"^  .  (136) 

r  r 

The  practical  significance  of  these  limitations  must  be  determined  experimen¬ 
tally  to  determine  a  meaningful  boundary  between  low- and  high-aspect  ratio. 
Some  experimental  correlation  has  been  obtained  on  a  very  low-aspect  ratio 
rectangular  wing  by  Fralich  and  Hedgepeth. 
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3. 


Supersonic  I^ifling  Surface;  Th(*ory 


a.  Int  roduclion 

For  many  applications  involving  low-aspect  ratio  planforuis, 

the  low  supersonic  regime  with  sul)sonic  leading  edges  is  next  in  imjjortance  after 

the  high  sulisonic  and  transonic  rc^giiTie.  For  the  purposes  of  achieving  a  unified 

format  for  calculating  the  AIC  matrices  by  the  relations  of  Eqs.  (35)  and  (  3f)), 

we  choose  the  supersonic  kernel  function  i)rocedure,  based  on  acceleration  po- 

52 

tential  as  developed  by  Cunningham,  Watkins,  and  Woolston  based  on  earlier 

5  3 

work  of  Watkins  and  Berman.  (See  Ref.  38  for  a  thorough  treatment  of  the 
various  methods  of  solving  supersonic  lifting  surface  problems.)  We  shall  now 
review  some  work  in  possible  alternative  methods  and  indicate  what  has  been 
done  to  develop  the  AIC  matrices  and  how  this  development  is  influenced  by  the 
type  of  a(irodynamic  solution. 


We  classify  the  solutions  of  interest  in  aer oelasticity  as  follows 

Class  1:  Closed  form  velocity  potential  solutions  for  specific  analytic 

forms  of  downwash  variation; 

Class  2:  "Box"  method  numerical  solutions  for  velocity  potential 

distribution; 

Class  3:  Pressure  loading  or  acceleration  potential  kernel  function 

solutions. 


b.  Discussion  of  the  Methods 

In  Classes  1  and  2,  the  integral  eejuation  has  been  inverted  ex¬ 
plicitly  to  ol^tain  the  velocity  potential  in  terms  of  the  downwash.  It  has  been 
shown  (e.  g.  ,  Refs.  34  and  38)  that  for  the  harmonic  case  [w  =  Re(-wexp(ikt)  )] 


0(x,  y ,  0:k) 


TT 


exp 

ikM^  (x-  4 ) 

7 

J 

®M.  f.  c. 

P  J 

cos  V(x  -  4)^  -  P^(y  -  n)^ 
V{x-e)^  -  p^(y-n)^ 


ds(e.n) 


(137) 
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For  the  so-called  "simple  planform"  (e.  g.  •  with  supersonic  edges)  the  solution 
is  complete  because  the  Mach  forecone  area  of  (x,  y),  ^  ^  ,  lies  completely 

on  the  planform.  For  subsonic  edges,  the  problem  is  complicated  by  the  fact 
that  is  unknpwn  in  the  region  between  leading  edges  and  the  Mach  cone 

emanating  from  the  wing  apex  (see  Fig.  8).  This  region  is  called  the  "diaphragm" 
region  for  which  the  boundary  Condition  is  the  vanishing  of  the  potential.  The 
box  method  solutions  of  Class  2  proceed  with  Eq.  (137)  formally,  and  compute 


the  diaphragm  downwashes  at  a  discrete  set  of  points  from  the  vanishing  of  the 
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potential  condition;  the  reader  is  referred  to  Pines,  Dugundji,  and  Neuringer, 

Zartarian  and  Hsu,^^  Zartarian,^^  and  Li^^  for  details.  Further  developments 

5  8 

have  been  made  for  computation  by  Ashley  and  Andrew.  The  numerous  static 
versions  (see  Etkin  and  Woodward,  for  example)  will  not  be  discussed  spe¬ 


cifically,  but  have  been  useful  in  problr*ms  of  optimum  camber  distribution.  Wc 


consider  the  prospects  for  box  solutions  in  the  collocation  AIC  format  briefly. 


Fig.  8.  Supersonic  Planform  Geometry 
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To  begin  with,  an  essential  point  with  the  collocation  format  is  the  lump¬ 
ing  of  loads  to  obtain  the  for ce -p r essur e  relationship.  The  box  method  assumes 
constant  pressure  and  downwash  on  each  "box"  so  that  such  an  element  is 
necessarily  a  much  smaller  subdivision  of  the  planform  than  is  the  AS.  associ- 
ated  with  collocation  control  point  force,  wherein  we  are  seeking  the  coarsest 
possible  grid  consistent  with  a  .good  representation  of  natural  vibration  modes 
shapes  and  frequencies.  Thus,  unless  chordwise  deformation  is  a  crucial 

fe.iture  of  the  aeroelastic  problem  in  question,  we  shall  seek  the  smallest  pos¬ 
sible  number  of  aerodynamic  control  points.  We  are  thur  led  to  the  notion 
that  a  ratiier  large  number  of  boxes  must  be  subassociated  with  a  collocation 
point  such  that 


F. 

J 


r. 

J 


AS. 


Ap 


(138) 


Now  Eq.  (138)  presents  no  difficulties  for  the  fixed  box  grid;  however,  the  con¬ 
sensus  experience  with  the  box  method  indicates  the  superiority  of  the  Mach 
box  approach,  wherein  the  spanwise  dimension  of  the  box  is  b^/P  and  the 

chordwise  dimension  is  b^  This  gives  the  same  number  of  boxes  covering  the 
entire  diaphragm  region  for  all  Mach  numl)ers  and  is  computationally  simpler 
to  handle.  On  the  other  hand,  as  M->l,P->0,  the  lateral  box  dimension  be¬ 
comes  proh ibitivedy  large  unless  tlie  basic  grid  is  proportionately  finer  (in 
which  case  the  computing  time  becomes  excessive).  Notwithstanding  this,  the 
solution  has  become  a  practical  one  as  generally  applied  to  the  computation  of 
modal  generalized  aerodynamic  forces,  for  say  M>1.  3.  Now  regarding  the 
collocation  method,  the  variable  location  of  Mach  box  grid  presents  a  great 
difficulty  in  the  logic  of  assigning  boxes  to  a  fixed  grid  control  point  so  that  we 
may  apply  Eq.  (1  38).  Furthermore,  to  satisfy  the  boundary  conditions  of  the 
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aerodynamic  problem  (especially  Ap  =  0  on  R),  the  kernel  function  approach  is 
superior.  In  view  of  the  importance  of  the  subsonic  edge  regime,  the  de¬ 
terioration  of  mechanical  accuracy  of  the  box  method  as  M->  1,  and  the  "float¬ 
ing  grid"  logic  problem,  we  turn  our  attention  to  the  methods  of  Class  3. 

It  is  desirable  to  have  a  unified  subsonic  -  transonic  -supersonic  lifting 

surface  approach  to  the  development  of  AICs.  For  this  purpose,  the  pressure 

kernel  function  approach  (Class  3  solution)  enables  use  of  much  of  the  computer 

coding  logic  developed  for  the  subsonic  program.  In  particular,  the  [W . .)  is 

i  i 

tlie  same,  as  is  the  relation  of  [  A  )  to  !D  1  [Eq.  (47b)l  .  As  will  be  seen, 
only  the  kernel,  the  loading  functions,  and  the  domain  of  integration  (S^^  ^  ^ 
instead  of  the  entire  planform)  are  changed;  hence,  the  ^^nm^  ma¬ 

trices  are  the  only  ones  changed.  Furthermore,  if  the  moments  of  the  chordwise 
pressure  loadings  used  in  the  [  development  employing  the  constrained 

minimum  vector  length  solution  [See  Eqs,  (71)  and  (7Z)]  are  developed  by  a 
numerical  quadrature  technique,  then  the  sarne[B^^]  matrix  program  may  be 
retained.'  For  a  large  number  of  chordwise  loads,  we  require  additional  con¬ 
ditions  of  the  form 


N.(y) 


V  B-*  h(x.) 

Z_  nm  '  j 


y+Ay/2 

dy 

y-Ay/2 


h(?.y)de 


(139) 


e . 


where  tlie  h(x)  represent  deflection  functions  whoso  generalized  forces  are  to  l)e 

preserved  by  the  equivalent  set  of  concentrated  loads  in  the  AIC  representation. 
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Ihe  alcove  approach  is  the  essence  of  the  current  line  of  development  of  AICs. 

The  remaining  Class  1  or  closed  form  solutions,  will  now  be  discussed. 

6)0  38 

In  this  category  are  the  works  of  Watkins  and  Berman  and  Miles  which  are 
based  on  low  reduced  frequency  solutions.  Ref.  60  considers  a  cubic  frequen¬ 
cy  expansion  for  a  wing  whose  deformation  is  of  the  form 


-'These  functions  may  be  any  functions  representative  of  the  type  of  mode 
shapes  thought  to  be  important  in  the  intended  application.  We  thus  find  some 
kinship  to  classical  modal  aeroelastic  analysis  here'. 


4a 


2  2 


h(x,  y'>  ~  Y.  Y. 


d  X 


(140) 


r=0  8=0 


Miles  (Ref.  38,  Sections  4.4  and  4.5)  illustrates  a  procedure  whereby  the  un¬ 
steady  solution  is  expressed  as  a  power  series  in  K  =  kM/O.If,  following  Miles, 


0(x/3,  y,  z,  0t)  =  0(x/e,y,z)e 


ikMt 


Hla) 


—  2  2 
0  =  §(x/0 ,  y ,  z)exp(ikM  x/0  ) 


;i41b) 


wherein  the  coordinate  and  time  are  made  dimensionless  by  reference  tojC  and 
i/V,  respectively. 

Then,  a  solution  is  found  of  the  form 


n=0 


For  the  low  frequency  case  the  solution  of  #  is  given  in  terms  of  the  steady 
state  term  (See  Ref.  38,  Eq.  4.5.2).  The  procedure  can  be  continued  for  the 
higher  As  a  starting  point,  the  solutions  of  Lance^^  provide  steady  state 
solutions  for  a  general  cubic  [^r  =  0(1)3,  s  =  0(  l)3j  in  Eq.  (140)  above.  Rodden 
and  Revell^^  have  developed  AIC s  for  narrow  delta  wings  utilizing  these  results. 
The  course  of  this  development  follows  closely  the  slender-wing  development 
described  in  Section  II.  B.  2,  and  the  loading  functions  can  be  expressed  as  a 
linear  combination  of  the  loadings  of  the  type  found  in  slender-wing  theory. 

The  coefficients  are  Mach  number  dependent,  involving  elliptic  integral  func¬ 
tions  . 

The  solutions  of  the  Class  1  type  have  the  difficulty  that  the  deflection 

interpolation  in  terms  of  Eq.  (  132)  "in  the  large"  can  be  inadequate  for  low- 

aspect  ratio  wings  with  extensive  chordwise  deformation.  However,  the  result 

of  Eq.  (  132)  can  be  improved  if  N.,  the  number  of  deflections,  exceeds  greatly 

J 

the  number  of  undetermined  coefficients  in  Eq.  (  140). 
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Before  concluding,  we  shall  mention  some  alternative  developments  of  the 
supersonic  kernel  function  procedure.  Ferguson  and  Guderly^^  propose  a 
solution  by  deriving  a  Volterra  integral  equation  of  the  second  kind,  whose 
domain  of  integration  is  confined  to  the  planform  and,  as  such,  provides  an 

9 

extension  of  Evvard's  method  to  unsteady  flow  for  the  case  of  subsonic  side 

edges  and  a  supersonic  trailing  edge;  however,  the  iterative  procedure  proposed 

to  solve  the  integral  equation  is  not  considered  as  simple  as  the  direct  matrix 

solution  of  the  Watkins’  acceleration  potential  kernel  function  method.  (For 

a  review  of  earlier  difficulties  in  attempting  to  apply  Evvard's  area  cancelling 

techniques  to  unsteady  flow,  see  ReTf.  38,  pp.  71-72.)  We  note  finally  that 

\/,  >  6*4 

Stanisic,  in  an  apparently  independent  series  of  investigations,  has  advocated 
the  use  of  the  acceleration  potential-kernel  function  approach. 

To  summarize  briefly,  the  supersonic  acceleration  potential  method  is 
considered  to  be  the  most  efficient  for  the  development  of  AICs  as  defined 
here  for  the  following  reasons: 

(1)  It  automatically  satisfies  the  pressure  boundary  con¬ 
ditions  in  the  diaphragm  region  R  (see  Fig.  8); 

(2)  It  utilizes  much  of  the  same  digital  computing  logic 
as  does  the  subsonic  problem  providing  a  more  nearly  unified  basis  for  lifting 
surface  theory  over  the  range  of  Mach  numbers; 

(3)  It  is  versatile  for  handling  arbitrary  planforms  and 
interference  effects  over  the  entire  range  of  reduced  frequency;  hence,  it  is 
suitable  for  practical  engineering  use  where  rapid  assessment  of  a  variety  of 
aerodynamic  planforms  is  an  important  capability  in  the  earlier  phases  of 
design; 

(4)  It  is  not  enslaved  to  a  particular  form  of  surface 
deformation  function  (as  with  the  Class  1  solutions); 

(5)  It  docs  not  suffer  additional  deterioration  of  accu- 
racy(beyond  the  implicit  limits  of  linear  potential  theory)  which  are  found  in 
the  Class  2  solutions  as  M-^1; 
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(6)  It  is  satisfactory  for  the  supersonic  edge  case 

(notwithstanding  the  possibility  of  an  explicit  solution  in  this  case);  hence,  the 
method  can  serve  as  a  logical  basis  for  a  unified  supersonic  lifting  surface 
theory  for  development  of  AICs.  (Further  results  of  Ref.  59  will  be  published 
subsequently. ) 

4.  Hypersonic  Lifting  Surface  Theory 

2c.  Introduction 

In  the  hypersonic  Mach  number  regime,  our  first  consider¬ 
ation  will  be  tlio  development  of  AICs  based  on  piston  theory.  The  later  con¬ 
siderations  will  discuss  the  various  classes  of  approximations  whose  validity 
depends  upon  the  hypersonic  similarity  parameter  where  is  the  effective 

local  body  surface  inclination  to  the  free  stream  with  due  account  of  angle  of 
attack  and  viscous  displacement  effects. 

b.  AICs  by  Piston  Theory 

The  piston  theory  has  been  mentioned  as  a  strip  theory  in 
Section  II.  A.  However,  as  a  point  pressure-downwash  relationship,  it  is 
more  appropriately  considered  as  a  lifting  surface  theory.  As  such,  the 
pressure  distribution  can  be  found  for  any  arbitrary  surface  deflection  mode, 
from  which  the  AICs  can  be  derived.  However,  this  solution  has  not  been  ob¬ 
tained  as  yet  in  general  form,  and  results  will  be  shown  here  only  for  the 
case  of  a  rigid  chord  with  a  control  surface.  The  derivation  of  the  oscillatory 
aerodynamic  coefficients  and  the  AICs  is  given  in  Ref.  65  and  will  not  be  re¬ 
produced  here.  It  differs  only  slightly  from  the  treatment  of  Ashley  and 
Zartarian,^^  in  that  the  third-order  pressure  coefficient  is  generalized  to 
account  for  sweep  and  steady  angle  of  attack,  and,  following  a  suggestion  of 
Morgan,  Huckel,  and  Runyan,  to  conform  with  the  second-order  quasi¬ 
steady  theory  of  Van  Dykc.^^  Therefore,  the  generalized  pressure  equation  is 
written  in  the  form 

p  -  Pq  =  C,(w/ap)  +  +  Cjiw/ap)^  .  (143) 

where,  for  the  usual  piston  theory. 
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=  -iKj/k  (146a) 

\o  "  -  -  ^^2/^  (146b) 

^pO  =  - -  '<^5  ■  (‘46c) 

~  -11^2/'^  (146d) 

'^aO  "  (146e) 

Mpp  =  -Kg/k^  -  i(K^  -  2K^?j^)/k  (146f) 

^hO  "  -*<^5  •  ^^4''h>/'^  (146g) 

■^aO  "  •  '<^6  ■  2K5?^)/k  (146h) 


^BO  =  -‘^5  -  2*^4 -  '(^6  ■  ^  4K^?^^)/k  (146i) 

where 

=  Xj^Zb  (117) 


K 

K 

K 

K 


1 


Z 


3 


4 


(1/m)|cj  -f  +  3C3M^(I^  -f 

(1/m)|cj  +  4C2MI2  +  3C3M^(2I^  +  a 
(4/3M)|cj  +  6C2MI3  +  3C3M^(3I^  + 

(1/m)  |Cj(1  -  ?|^)  +  2C2MJJ  +  3CjM 


a  ^1 
0 


(  148a) 

(  148b) 
(  148c) 

(  148d) 


47 


1 


( 144a) 
( 144b) 


=  (Y  +  l)/4 
C3  =  (Y  +  1)/12 


and,  for  the  quasi- steady  theory, 


C 


1 


2  A 

sec  A 


(  144c) 


(  145a) 


C 


2 


M^(Y  + 


2  2  2 
1)  -  4  sec  A  (M  -  sec  A 

Z 


4(M 


-  sec 


'A 


(  145b) 


Cj  =  (Y  +  1)/12 


(  145c) 


where  the  value  of  is  taken  from  piston  tlieory  since  Van  Dyke  gives  only  a 

^  25 

second-order  solution.  Ashley  and  Zartarian  have  shown  that  the  piston 

2  2 

theory  is  applicable  if  any  of  the  conditions  M  ;:^1,  Mk»l,  or  k  »1  are  met. 
We  see  that  if  the  reduced  frequency  is  low,  the  Mach  number  must  neces¬ 
sarily  be  high;  liowever,  if  the  reduced  frequency  is  large,  the  Mach  number 
need  not  necessarily  be  large,  and,  in  fact,  it  could  be  transonic  or  even  sub¬ 
sonic.  From  this  it  is  apparent  that  any  correction  introduced  to  bring  piston 
theory  into  line  with  linearized  supersonic  theory  must  be  considered  as  a  low 
frequency  approximation.  The  use  of  Eq.  (145)  extends  the  lower  l  lach  num¬ 
ber  limit  of  piston  theory  at  low  frequencies.  Eq.  (145)  is  seen  to  be  the 
more  general  result,  since  if  sec  A  is  taken  to  be  zero,  Eqs.  (  144)  result,  and 
no  correction  is  made. 


The  derivation  of  the  oscillatory  aerodynamic  coefficients  referred  to 
the  leading  edge  (i.e.  ,  both  motion  and  pitching  moment)  leads  to  the  following 
value  s: 


=  -iKj/k  (146a) 

^qO  =  -  (146b) 

^-pO  "  -  1(1^5  -  2K4^h>/'^ 

^hO  "  -IK2/*"  (1‘16'1) 

^ao  "  ('“l^e) 

MpQ  =  -K^/k^  -  i{K^  -  2K3?^)/k  (146£) 

\o  =  -^<^5  -  ^^4'’h>/’'  (146g) 

■^aO  "  -<^^5  -2K4?h)/'^^  -  i(1^6  '  (146h) 


Tgo  =  -(K5  -  2K^5,,)/k^  -  i(K^  -  4K3?^  +  4K4^j^^)/k  (1461) 

wh  ere 

=  Xj^Zb  (117) 


K 

K 

K 

K 


1 


Z 

3 


4 


(1/n4)|cj  + 

(1/m)|cj  +  4C2MI2  +  3C3M^(2I^  +  a 
(4/3M)|cj  +  ^>C2Ml3  +  3C3M^(3I^  + 
(1/m)  |Cj(1  -  ^1^)  +  +  3C^M 


-  ^h) 


(  148a) 

(  148b) 
( 148c) 

(  148d) 
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Kg  =  (1/m)|Cj(1  - 


^6  = 


r  ,  2 

zl 

2J.  +  a  (1  - 

L  5  o' 

h  iJ 

- 

r  ,  2 

3  ~ 

|_3J6  +  (1  - 

148f) 


and  the  thickness  integrals  for  the  airfoil  having  the  symmetrical  thickness 
distribution  2g(x)  are 


/2b  ,  , 

(l/2b)^  “  1»2,3 

{l/2b)*^  ^(^)^dx,  for  n  =  4,5,6 

0 

{l/2b)*^  j  x^  for  n  =  1,2,3 


l/2b)^ 


dx,  for  n  =  4,5,6 


The  thickness  integrals  are  evaluated  in  Ref.  65  for  an  airfoil  cross  section 
idealized  into  a  parabola  from  the  (sharp)  leading  edge  to  the  point  of  maximum 
thickness,  another  parabola  to  the  control  surface  hinge  line,  and  a  straight 
line  to  the  (possibly  blunt)  trailing  edge. 

The  development  of  the  AICs  in  Section  II.  A  assumed  the  oscillatory 
coefficients  were  referred  to  the  quarter- chord.  Since  the  present  coefficients 
are  referred  to  the  leading  edge,  an  alternative  derivation  of  the  AICs  for  the 
strip  is  necessary.  This  is  also  given  in  Ref.  65,  assuming  again  the  control 
points  to  be  at  the  quarter-chord,  hinge  line,  and  trailing  edge,  and  the  re¬ 
sulting  expression  is 
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4(b/b^)^(£,y/s) 


(1  +  b/2d)  -b/d 
-b/2d  b/d 
0  0 


(b/c  K3b/2d-l) 

<x 

-(b/cj(3b/2d) 

b/c 

a 


X 


E  ^ 

L 

hO 

aO 

0  0 

'^hO 

aO 

^60 

^hO 

H 

P 

o 

^BO 

(1  +  b/2d) -b/2d  0 

-b/d  b/d  0 

b/d  -(b/d+b/c  )  b/c 
a  a 


(151) 


which  in  the  absence  of  a  control  surface  reduces  to 


Cj1=  4(b/b^)^Ay/s) 


(1  +  b/2d)  -b/d 

^hO  '"aO 

-b/2d  b/d 

.^hO  ^ao. 

(l+b/2d)  -b/2d" 
-b/d  b/d 


.(152) 


We  have  illustrated  a  technique  to  extend  the  lower  Mach  number  limit  of 

piston  theory.  The  question  remains  as  to  the  upper  Mach  number  limit.  This 

question  has  been  discussed  by  Morgan,  Runyan,  and  Huckel,^^  and  by  Miles. 

69 

Experimental  correlation  has  been  obtained  up  to  M  =  6.  86.  Miles  argues 
that  the  basic  philosophy  of  piston  theory  is  valid  if  changes  in  entropy  associ¬ 
ated  with  tlie  unsteady  flow  can  be  neglected.  Then  the  perturbation  pressure 
p'  may  be  calculated  from  the  [da  no -wave  relation 


p'  =  paw'  (153) 

where  p  and  a  denote  the  local  density  and  sonic  speed  for  the  steady  flow, 
and  w'  is  the  unsteady  component  of  the  downwash  velocity.  It  is  emphasized 
that  p  and  a  need  not  be  based  on  the  assumption  of  an  ideal  fluid.  In  short, 
Eq.  (153)  states  that  the  unsteady  solution  is  known  when  the  steady  one  is  ob¬ 
tained,  provided  the  changes  in  entropy  due  to  unsteadiness  may  be  neglected. 
Miles  suggests  procedures  for  calculating  p  and  a  when  the  usual  piston 

theory  li.  e.,  Eq.  (143)  becomes  inadequate.  Morgan,  Runyan,  and 
68^ 

Huckel  suggest  that  a  probable  upper  limit  of  validity  of  piston  theory  is 
about  Macli  15.  Beyond  this,  the  problems  of  entropy  changes,  di sso- jiation, 
and  boundary  layer  interaction  must  be  considered. 
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c.  General  Hypersonic  Inviscid  Considerations 

In  the  liypersonic  regime,  a  given  blunt-nosed  slender  body 
or  wing  can  be  thought  of  as  possessing  several  distinct  local  flow  regimes  de¬ 
pending  on  the  local  hype rsonic  similarity  parameters  (see  Ref.  70,  Sec.  1.  3). 

(154) 


K  =  Me^(x) 

b 


=  Ma(t) 


155) 


where  is  the  local  streamwise  slope  of  the  body.  If  we  specialize  to  the 

case  of  small  unsteady  motion,  we  may  write  for  longitudinal  motion 


a  =  -t  h(x,t)/v  +  a(x-XQ) 


where 


h/V 


,  ai/v 


«  1 


156) 


(157) 


Then  in  this  case  rw  and  becomes  fixed  for  a  given  steady  state  vehicle 

a  ttQ 

attitude  a  .  Now  for  any  blunt-nosed  vehicle  0|^  varies  from  90  degrees  to 
nearly  zero.  We  must  also  add  d  ^'/dx  »  the  boundary  layer  displacement 

thickness  gradient,  to  the  slope,  and  thus  the  effective  local  streamwise 
slope,  9^  +  d6  /dx,  is  always  finite. 

We  now  consider  a  typical  hypersonic  wing  section  as  shown  in  Figure  9. 
For  convenience,  we  note  that  Hayes  and  Probstein  classify  the  hypersonic 
regimes  as  follows: 

Thin  shocklayer  or  ’’strong  shock”  (nose)  regime 


(1) 


MP,  (x)>>  1 

b 


(158) 


(Z)  Small  density  ratio  regime 

=  "2n^"on> 
e«l 


159) 
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SECTION  B-B 


Fig.  9.  Hypersonic  Planform  and  Airfoil  Section 


(3)  Slender  body  regime 

e^(x)«l 


(160) 


(4)  Linear  regime 


Me^(x)«l 


(161) 


In  the  nose  regime,  the  simplest  approach  for  aeroelastic  purposes  is  the 

Newtonian  impact  theory.  Applications  of  Newtonian  theory  to  unsteady  flow 

have  been  given  by  Hayes  and  Probstein,^^  Kennet,^^  and  Tobak  and  Wehrend.^^ 
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A  typical  steady  flow  application  is  given  by  Jackson.  Not  all  thin  shock 

71 

layer  problems  are  adequately  described  by  Newtonian  laws.  Kennet  has 
studied  unsteady  perturbations  in  the  stagnation  point  neighborhood  within  the 
shock  layer  due  to  plunging  oscillations,  and  shows  some  severe  departures 
from  quasi-steady  Newtonian  theory  at  M  =  5.  8.  In  this  approach,  he  derives 
a  third -order  linear  partial  differential  equation  for  u,the  unsteady  meridian 
tangential  velocity  component  in  the  shock  layer,  (by  manipulation  of  the  mo¬ 
mentum  and  continuity  equations),  which  he  solves  by  a  power  series  in 


r|  -  y/6 


(see  Figure  10).  He  finds  shock  oscillation  amplitudes  related 
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SECTION  A-A 


A 


Fig.  10.  Blunt  Body  Shock  Layer  Geometry 


to  the  plunging  motion  and  surface  pressure  derivatives  with  respect  to  h,  h, 

04 

and  h.  This  method  requires  knowledge  of  the  steady  flow  within  the  shock 
layer,  which  is  rotational  and  is  nonuniform,  and  requires  satisfaction  of  the 
Rankine -Hugoniot  relations  at  the  bow  shock.  Tliese  features  are  similar  to 
those  discussed  with  respect  to  the  transonic  problem;  however,  here  the 
vorticity  is  of  tlie  essence  in  describing  the  flow  field  inside  the  shock  layer. 
From  the  point  of  view  of  tlie  aeroela  stician,  interest  centers  on  the  trans¬ 
verse  loading  and  therefore  on  the  loading  on  the  part  of  the  body  where 

G,  <<1,  but  where  still  MG,  is  large  for  high  hypersonic  speeds.  For  this 
b  o  ^ 

situation,  Newtonian  theory  should  be  useful,  and  the  stagnation  point  diffi¬ 
culties  present  no  practical  problem  since  these  pressure  fluctuations  are 
directed  chordwise,  against  a  rather  solid  nose. 

d.  Newtonian  Theory 

The  Newtonian  pressure  coefficient  is  given  by 
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Fig.  11,  Orientation  of  the  Unit  Normal  Vector 


Then  for  two-dimensional  flow,  the  lower  and  upper  surface  pressure  coeffi¬ 
cients  are  given  bv 


Cp  ^  -  2  sin^P^i  (Wj^/V)  sin  26^  +  cos^P^  (165a) 


For  bodies  of  revolution  we  have 


C  “2  sin^0  +  (\^,  /V)  sin  2B  sin  9  +  (w,  /V)^  cos^  0.  sin^  0  (165b) 

nL.bb  b  b  bj. 


We  have  the  requirement  of  positive  pressure  for  all  Newtonian  theories 

C  0 
P"" 

For  practical  purposes,  Newtonian  theory  can  be  improved  by  using  Lees' 
empirical  suggestion  and  replacing  the  factor  2  by  the  stagnation  point  pres¬ 
sure  coefficient  as  is  done  in  steady  flow.  This  empirical  form  is  simpler  to 
use  than  the  more  logically  appealing  Busemann  correction  for  centrifugal 
force  which  unfortunately  can  lead  to  negative  pressures  (see  Ref.  70,  Ch.  3). 

e.  Recent  Methods 

An  alternative  thin  shock  layer  approach  has  its  roots  in 
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"blast  wave  theory.  "  Recently  Zartarian,  Hsu,  and  Ashley  have  developed 

a  variational  approach  for  the  equivalent,  two-dimensional  piston  problem  in 

the  plane  transverse  to  the  body  axis.  This  theory  is  aimed  at  the  neighbor - 

hood  MP,  ~  1  to  supplement  their  extension  of  tlie  well-known  shock  expansion 

theory  to  unsteady  flow  based  on  the  quasi-steady  assumption.  Because  many 

liypersonic,  unsteady  aeroelastic  plienomena  will  occur  at  very  low  reduced 

frequencies,  the  quasi-static  analyses  appear  more  attractive  than  at  transonic 

speeds.  In  this  category,  tlie  sliock  expansion  method  as  extended  by  Zartarian, 

H.su,  and  Ashley  is  especially  pertinent,  for  the  effects  of  the  nose  shock  can 

be  handled  quite  easily.  These  investigators  find  distinct  improvements  over 

the  widely  used  piston  theory  for  the  range  ^  ^  ^  although  piston  theory 

is  fairly  adequate  up  to  M  9,  —  0.  7. 

b 

The  above-mentioned  theories  all  have  the  simplicity  of  a  point  relation 
between  pressure  and  downwash,  and  offer  few  difficulties  to  aeroelastic 
applications.  All  of  the  theories  show  the  importance  of  thickness  effects 


(166) 
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which,  as  previously  mentioned,  should  include  the  viscous  displacement  thick¬ 
ness  We  next  consider  the  viscous  interaction  problem. 

f.  Viscous  Interactions 


The  prol)lem  of  viscous  interactions  has  received  great  atten¬ 
tion  in  steady  hypersonic  flow  theory,  since  aerodynamic  control  considerations 
are  concerned  with  accurate  predictions  for  Mach  numbers  and  altitudes  much 
higher  than  those  at  which  maximum  dynamic  pressures  are  expected.  In  these 
cases  of  low  Reynolds  numbers  the  problems  are  very  serious,  and  we  find, 
for  example,  that  the  concept  of  a  sharp  leading  edge  is  completely  illusory; 
we  have  at  least  a  "half-power"  body  since  6 'I'/x ^ ( y- 1)M  /  yRe^.  Now 
at  verv  high  Mach  numbers  and  (6„,)  V vd,  JN  ,  and 

the  shock  layer  and  boundary  layer  tend  to  merge  at  low  Reynolds'  numbers. 


However,  if  the  leading  edge  diameter  Reynolds  number  Re^  is  sufficiently 
large,  then  the  boundary  layer  thickness  6  is  much  less  than  the  shock  layer 


thickness  ^32^' 
According  to  Lees, 


and  the  inviscid  flow  dominates  the  stagnation  region, 
this  condition  is  met  if 


Under  such  circumstances,  the  analysis  of  the  type  conducted  by  Kennet  is 
adequate  for  the  stagnation  point  vicinity. 

We  have  already  seen  the  importance  of  thickness  on  the  pressure  distri¬ 
bution  in  the  treatment  of  piston  theory;  we  consider  therefore,  following 
Z.artarian,  Hsu,  and  Ashley,  that  to  the  first  approximation  the  addition  of 
the  displacement  thickness  6'''  accounts  sufficiently  for  the  effects  of  viscosity. 
Tiiis  factor  becomes  important  at  very  high  Mach  numbers  since  the  laminar 
boundary  layer  contribution  to  is 


A(M9^^) 


(Y-l)M^ 


X 


(168) 
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where  x  is  the  viscous  interaction  parameter.  It  is  well  to  recognize  that 
the  calculation  of  displacement  thickness  at  high  Mach  number  requires  some 
familiarity  with  the  field  of  aerothermodynaniics.  We  therefore  illustrate  the 
necessary  formula s .with  example  calculations.  The  aeroelastician  is  primar¬ 
ily  concerned  with  high  dynamic  pressure,  which  keeps  the  considerations  well 

4 

within  tlie  valid  regime  of  continuum  boundary  layer  tlieory  (say  >  10  , 

see  Ref.  70,  Ch.  10,  or  Ref.  40,  Ch.  14).  L.et  us  express  Re^  in  terms  of 
dynamic  pressure  and  free  stream  Mach  number. 


RGx  =  (169a) 

=  (l/^)pQV^(2x/M<y^^Up)  (169b) 

=  q(2x/MaQUQ)  (169c) 

where  tlie  zero  subscript  refers  to  the  free  stream  condition.  For  instance,  if 

-7  3  3 

Up,  =  4  X  10  Ib-sec/sq  ft,  a^^  =  10  fps,  and  q  =  10  psf,  then  Re  = 

5x10  (x/M),  which  for  M  =  10,  x  =  1  ft,  yields  Rc^  =  5  x  10  .  Since 

transition  Reynolds'  numbers  are  of  this  order  (see  Ref.  77,  p.  438),  much 

of  the  boundary  layer  will  be  turbulent. 

For  a  flat  plate,  the  displacement  thickness  is  given  by  (Ref.  77,  p.  537) 


=  (^'Vp)(C^/C^.)(C^.x/2)  (170) 

wliere  the  momentum  thickness  9  is 

P  =  (Cj,/Cp..){Cj,jx/2)  .  (171) 

As  an  example,  we  consider  values  for  C„/C„.  and  C  based  on  experimental 

/  Jt*  Jt*  1  r  1 

results.  Taking  Re  =  10  and  C  =  0.  003  (see  Ref.  77,  p.  538),  and  using 

^^78 

the  v'-alucs  of  H  =  given  by  Tucker  for  a  one-seventh  power  boundary 

layer  v'elocity  profile,  we  construct  Table  1  for  an  insulated  plate. 
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Table  1 


M 

H  = 

Cp/Cp.  (Ref.  40) 

6  /x«d6  /dx 

0.  1 

1.  29 

1.  00 

0.  00193 

1.0 

1.  73 

0.  93 

0.  00242 

5.  0 

12.  31 

0.  37 

0.  00691 

10.  0 

44.  92 

0.  16 

0.  0108 

We 

see  from  Table  1  that 

the  turbulent  displacement 

effect  at  M  =  10  : 

five  times  greater  than  M  =  0.  1.  For  actual  bodies,  we  must  rely  on  the 
aerothermodynamicist  to  compute  the  proper  values  of  6  considering  the 
presence  of  pressure  gradient,  surface  roughness,  heat  transfer,  ablation, 
real  gas  thermodynamics,  body  geometry,  and  any  other  practical  consider¬ 
ations.  Such  calculations,  however,  are  required  for  any  design  in  order  to 
assess  the  heat  transfer  problem.  In  the  laminar  boundary  (roughly 

Re  <  5  X  10^  or  Rep  <  480)  simple  estimates  are  often  possible  using  the 

^  ^  79 

local  similarity  concept  introduced  by  Lees  (see  Ref.  80  for  a  recent  review 

of  this  method).  However,  in  the  laminar  regime,  the  leading  edge  interaction 

problem  requires  a  successive  approximation  solution  to  the  boundary  layer 

growth  calculation  (see  Ref.  70,  Ch.  9,  and  also  Ref.  81).  For  the  more 

practical  case  of  a  noninsulated  wall,  the  turbulent  skin  friction  coefficient 

can  be  calculated  roughly  according  to 

(17^) 


if  tlie  viscosity  is  proportional  to  T  '  We  might  take  as  the  average  of 

tlie  wall  temperature  and  the  maximum  temperature  as  a  reasonable  estimate, 
although  other  choices  have  been  made  (see  Ref.  70,  pp.  328-9).  From  the 
definitions  of  the  boundary  layer  displacement  and  momentum  thickness 

(l-nu/noV)d(y/6)  (173) 
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By  differentiation,  we  obtain 


^  =  0.8  0-074x~°-^  f i_  \ 

dx  (Re-)°-2  Cp..  J 

=  0.86* /X  .  (178b) 

sjs 

Thus  d<^  /dx  —  0.  016  for  the  noninsulated  smooth  plate  example.  On  a  rough 
plate  —  constant,  and  dfi  /dx  “si  iS  /x  .  We  thus  find  a  rather  substan¬ 

tial  displacement  effect,  equivalent  to  a  1,6  percent  thickness  ratio  for  the 

example  conditions,  Re  =  10^,  M  =  ZO,  T  /T^  =  6  for  a  turbulent 
^  X  w  0 

boundary  layer.  For  these  conditions,  the  turbulent  boundary  layer  inter- 

75 

action  parameter  is  MdA  /dx  =  0.  3Z  .  Zartarian,  Hsu,  and  Ashley  have 
noted  the  deterioration  of  piston  theory  for  >  0.  3  so  that  in  the  example 

shown,  the  viscous  effect  is  significant  in  that  the  limits  of  second-order 
piston  theory  are  approached. 

The  above  calculations  would  be  somewhat  pessimistic,  since  dissocia¬ 
tion  of  oxygen  begins  at  T  ~  3500°  K  and  the  energy  absorbed  would  probably 
limit  temperatures  to  this  vicinity.  Furthermore,  the  vibrational  degrees  of 
freedom  for  oxygen  and  nitrogen  become  appreciably  excited  at  lower  temper¬ 
atures  (approximately  ZZ50°K  and  3370°K,  respectively),  and  the  ratio  of 
specific  heats  approaclies  4/3  above  tlie  higher  temperature,  thus  somewhat 
reducing  the  temperature  because  of  the  locally  higher  heat  capacity  near 
T  ;  a  complete  real  gas  calculation  might  show  a  smaller  value  of  6  /R  . 

In  this  connection,  Whalen  has  shown  that  for  air,  nonequilibrium  thermo¬ 
dynamic  effects  become  important  at  speeds  above  10,  000  feet  per  second  for 
altitudes  above  150,000  feet,  and  that  above  ZOO,  000  feet  the  gas  composition 
along  any  streamline  can  he  considered  "frozen"  at  its  equilibrium  value  on 
tlie  downstream  side  of  the  sliock  (see  also  Ref.  83).  Whalen  shows  that  non- 
equilibrium  behavior  considerably  affects  the  blunt  nose  influence  on  the 
steady  plate  flow  downstream  (and  hence  the  local  thermodynamic  state),  but 
does  not  seriously  affect  the  leading  edge  viscous  interaction.  For  simple 
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estimates  of  the  combined  effects  of  viscosity  and  blunting,  see  Bertram  and 
8 1 

iBlackstock. 

Returning  to  the  aeroelastic  problem,  we  consider  the  dynamic  pressure 
in  relation  to  Re^,  nM,  and  atmospheric  properties.  For  q  =  1000  psf,  M  =  ZO, 

Re  =  10^  ,  tlien  from  Eq.  (169c),  we  have  x  =  4  feet;  this  condition  exists  at 

^  •  Z 

an  altitude  of  145,  000  feet  (since  q  =  y  n  M  /Z  »  we  find  p^  =  3.  57  psf).  A 

^0  ^84 

recent  treatment  of  super  circular  orbital  lifting  re-entry  vehicles  shows 
altitudes  as  low  as  100,  000  feet  at  V  =  30,  000  feet  per  second  near  the  pullout. 
In  sucli  a  regime,  we  would  reach  very  high  dynamic  pressures  at  very  high 
Mach  numbers,  and  the  viscous  displacement  would  be  extremely  important. 

In  summary,  we  see  that  thickness  effects  are  an  essential  feature  of  the 
hypersonic  airloads  problem,  and  these  must  include  viscous  corrections. 
Furthermore,  any  given  hypersonic  body  is  describable  as  having  several 
regimes  of  the  local  hypersonic  similarity  parameter,  ,  which  dictates 

the  choice  of  lifting  pressure  approximation  to  be  used  locally. 

C.  AICs  for  Slender  Bodies 

1 .  Slender -Body  Theory 

The  equations  given  by  Bisplinghoff,  Ashley,  and  Halfman  (Ref.  34, 
p.  418)  provide  a  convenient  basis  for  deriving  the  AICs  of  a  slender  body.  The 
vertical  force  acting  per  unit  length  of  the  body  is  the  reaction  to  the  substantial 
rate  of  change  of  the  momentum  of  the  virtual  mass  per  unit  length  of  the  body 


dx 


(179) 


where  the  momentum  of  the  virtual  mass  per  unit  length  is  found  from  the  local 
cross-sectional  area  and  the  downwash 


dx 


pSw 


=  pS 


(  180a) 
(  180b) 
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If  we  assume  harmonic  motion  and  substitute  Eq.  (180b)  into  Eq.  (179), 
we  obtain 


dx 


X 


+  iwh 


J 


-  itopS  V  +  iwhj 


(181) 


To  obtain  the  force  on  a  specified  length  of  the  body,  it  is  necessary  to  inte¬ 
grate  Eq.  (181)  over  that  length.  We  consider  the  body  to  be  divided  into  a 
number  of  sections,  not  necessarily  of  equal  length,  as  shown  in  Figure  12. 
For  the  n'th  section,  the  control  point  is  taken  at  the  midpoint  of  its  length 

1  =  X  -  A  /2  and  has  the  cross- 
n  n 


The  aft  end  of  the  section  is  located  at  x 


^  “  2 

section  area  S  1  ;  the  forward  end  of  the  section  is  located  at  x  ,1  = 

"  -  2  "  +  Z 

X  +  A  /2  and  has  the  cross-section  area  S  1  .  Carrying  out  the  integration 
n  n  n  -1- 

of  Eq.  (181)  for  the  n'th  section  yields  the  following 


X  1 

n  J  dx 

X  1 
n  — 

2 


182a) 


=  -  pV 


V  ^  +  iwh 
dx 


->x  ,  1 

"■"z 

X  1 
"-Z 


-  lo)p 


^n+-  /  \ 

r  2  s  V  4^  +  icoh  dx 
J  \  dx  / 

X  1 

"-2 


(182b) 
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X..  .  X. 


Fig.  12.  Slender-Body  Geometry  for  AICs 


F  = 
n 


pv 


+  iw 


S(Vh'  +  iuh)  dx 


(  182c) 
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We  resort  to  Lagrangian  interpolation  for  the  evaluation  of  the  terms  in  Eq. 
(182c).  For  the  first  two  terms  we  choose  parabolic  interpolation 


h  =  h 


(x  -  x^)(x  -  Xn^,) 


(x-x„.))(x-x^^l^ 

'  ’'n+l' 


(x  -  Xn.iXx  -  Xr) 

K+1  -  -  V 


For  evaluation  of  the  integral  we  use  only  linear  interpolation 


h  =  h  +  h’  (x  -  X  ) 
n  n  n 


(1/A)S  l(x,l-x+S,lx-x  1 
'  '  n'  n--~  n+^  n+-  n- - 

2  2  2  ,  Z 


The  three  point  interpolation  leads  to  three  elements  in  each  row  of  the 
AIC  matrix.  Therefore 


F  =  poj^b^s  (C,  ,  h  1  +  C,  h  +  C,  ii  ^  .li 

n  r  hn,n-l  n-1  hn,  n  n  hn,  n+1  n+T 
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By  substituting  Eq.  (186)  into  Eq.  (l8Zc),  after  evaluating  the  integral  of  Eq. 
(182c)  by  means  of  Eqs.  (184)  and  (185),  we  obtain  the  following  relation  for  the 
AICs 


C,  ,  h  ,  +  C,  h  +  C.  1  1  ^  0.1 
hn,n-l  n- 1  hn,  n  n  hn,n+l  n+1 


(1/kS)^ 


Is  .  1 

h  1  - 

S  1  h' 

1 

1  n+  — 

n+~ 

n - n- 

1  ' 

i  2 

2 

2 

2 

h  AV  +  (A^/12)  S  ^  1  -  S  1 
n  n  n  n+^  n  — 

•  \  2  2 


n 


-  i  (kj,/b^). 


S  ^  1  h  ^  1  -  S  1  h  1  +  h’  AV„ 

n+—  n+— '  n- —  n- —  n  n 

2  2  2  2 


(187) 


where  is  tlie  volume  of  the  n'th  section. 

If  we  evaluate  the  deflections  in  Eq.  (183)  at  stations  x  1  and  x  1  , 

and  then,  by  differentiating  Eq.  (183),  evaluate  the  slopes  at  stations 
X  1  ,  X  ,  X  J  ,  we  may  then  place  the  riglit-hand  side  of  Eq.  (187)  in 

n-’^-  nnh^' 

terms  of  the  control  point  deflections,  h  , ,  li  ,  and  h  ,  ,  .  If  we  carry  out 

n  - 1  n  n  + 1 

this  substitution,  the  AICs  are  found  by  identifying  the  coefficients  of  the  con¬ 
trol  point  deflections  on  both  sides  of  Eq.  (187).  The  following  results  are 
obtained 
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\ 


where 


(l/k^s 

r 


C 


hn,  n-  1 


^-l)  [-  -  ^„+l)  -  2S„.l  (3A_^  +  A^^j) 


_2  „  2. 


-  S  I  -  S  iUa  +A  ,) 

n  +  -  "'jl  "  f'+l 


+  i  (k  /b  ) 
r  r 


S  1  A  A^,+S  lA  (2  A  +  A 
n  +  -  n  n+1  n--  n  '  n  n-f  T 


4  2AVn  (A_^  +  An^l)|}  (188) 


D  =  (A  4  A  )(A  ,  4  2A  4  A  ,) 

n-1  n-1  n'  n- 1  n  n+1' 


(189) 


^hn,  n 


(^n-1  +  -  ^n4l)  -  ^S,,.  ^  (A^_  j  -  2A^  -  A^^j) 


2  „  2 


+  (k  /b  )  D  AV  +  (A  /6) 
r  r  n  n  n 


S  ^  1  -  S  1 1  (A  ,  -  A  ) 
n  +-  n  -  —  '  n+1  n-  1 ' 


-  (Vl  (^^n-1^4l) 

^  2 


4  2AV  (A  -  A  ,) 
n  1)41  n-1' 


J 


(.1  90) 


OS 


where 


(191) 


D 

n 


(A  ,  +  A  )(A  +A 
n- 1  n  n  n+1 


The  above  expressions  are  applicable  for  all  intermediate  sections  of  the  body, 
sections  which  are  centrally  located  as  far  as  the  interpolation  and  differentia¬ 
tion  are  concerned.  The  exceptions  arc  tlie  first  and  N'th  sections.  The 
counterpart  of  Eq.  (187)  for  the  first  section  is 


bb 


^hl,  1  1  *"hl.  E  •’Z  Si,  3  ’’3  ■ 


(l/k^s)  < 


2  ,.  2. 


-  S,  h'  +  S,  h'  +  (k  /b^) 
3  3  11  r  r 

2  I  22 


hj  AV,  +  (Aj/12)(S3  -  Sj)  h'j 

2  2 


-  i  (k^/b^XSj  hj  -  Sj  hj  +  hj  AVj) 
2  2  2  2 


(194) 


Carrying  out  the  evaluation  of  the  appropriate  deflections  and  slopes  from 
Eq.  (183)  in  terms  of  the  first  three  control  point  deflections  leads  to  the  fol¬ 
lowing  coefficients  for  the  first  row  of  the  AIC  matrix 


"hi.  1 


(1/k^s  Dj)  {3A^  +  A3)  -  2Sj  (4A^  +  3A2  +  A3) 


+  (k^/b^) 
r  r' 


Dj  AVj  -  (Aj/6)(S3  ■  Si)(2Aj  +  3A2  -  A3) 


2  2 


83  A^  (2A2  +  A3)  -  Sj  (2Aj  +  A2)(2Aj  +  2A2  +  A3) 

L  2  2 


-  2AVj  (2Aj  +  3A2  +  A3) 


(195) 


wh  ere 


D,  =  D  ,  with  n  =  Z 
1  n  - 1 
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where 


wliere 


"hi. 2 


(l/k^8  D^)<ZS^  (A^  -  ZA^  -  A3)  +  2Sj  (3Aj  +  ZA^  +  A3) 


+  (kJ/b^)(A^/6)(S3  -  Sj)(Aj  +  2A2  +  A3) 


z  z 


-  i  (k  /b  ) 
r  r 


S3  A^  (ZA^  +  A3)  +  Sj  Aj  (ZAj  +  ZA^  +  A3) 
+  ZAVj  (A^  +  ZA^  +  A3)|^ 


(196) 


D-,  =  D  with  n  =  Z 

Z  n 


"hi,  3 


r 


(l/kS  D3)<j-  2S3  (Aj  -  A^)  -  2Sj  (3Aj  +  A^) 
2  2 

-  {kJ/bJ)(Aj/6)(S3  -  Sj)(Aj  +  A^) 


+  i  (k  /b  ) 
r  r 


2  2 

SsAj  A^  +  SjA,  (2A,  +A3I 


+  2AVj  (Aj  +  A^) 


(197) 


D-,  =  D  .  with  n  =  Z  . 
3  n  +  1 
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Similarly,  the  counterpart  of  Eq.  (187)  for  the  N'th  section  is 


^liN,  N-2  *'n-2  ^hN',  N-1  ’'N-1  ^  Sn,  N 


('/V)  -SN4‘'N4+SN4Vi 

I  Z  z  z 

+  (k^/b^)  (h^  AV-,  +  (A^,/12){S„^  1  -  S„  1)  h' 

'  r  r'  N  N  '  N  ''  N+—  N —  N 

Z  Z 

-  i  (k  /b  i  1  -  S.,  1  h..  1  +  h’  .  „ 

r  r''  N+y  N+—  N- y  N- y  N  AV^, 

Cs  Ls  ^  M  IN 


(198) 


and  we  obtain  the  following  coefficients  for  the  last  (M’th)  row  of  the  AIC  matrix 


"hN,  N-Z 


Dn-2)|-  +  ^N-l) 

+  (kJ/bJxA^/bXS^^^-  Sj,  1)(A^  + A^_j) 

'  z  z 

-  i(kj./b^)  (2Aj^  +  Aj^  j)  +  1  A^^  Aj^_  j 

Z  —  z 

+  2AV^  (A^H-An.i)]. 


(199) 


where 


D-.  ,  =  D  ,  with  n  =  N-1 
N-Z  n-1 
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'"hN,  N-1  ' 

(l/kS  (3A^  +  2A^_j  +  A^_2)  +  2S^  1  (A^  -  2A^,i  -  A^,^) 

N  2  2 

-  (kJ/b2)(A^/6)(S^  1  -  S^  l)(Aj^  +  2A^,j  +  A^_2) 

+  i(^/b,)Js^,  jA^  (2A^  +  2A^_j  +  A^_2) 

■*■  ^N--^N  ^N-2^ 


+  2AVj^  {A^  +  +  Aj^_2)1 


(200) 


wh  ere 


^N-1  ~  ^n  ^  ~ 
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we 


Before  going  on  to  a  discussion  of  second-order  slender  body  theory, 
note  that,  according  to  Miles,  the  slender -body  theory  presented  above  may 
be  applied  to  finned  vehicles  if  an  effective  cross-section  area  in  the  finned 
region,  whose  geometry  is  shown  in  Fig.  13,  is  taken  as 


S 

e 


(203) 


Fig.  13.  Cross  Section  of  Finned 
Region  of  Slender  Body 


2.  Second-Order  Effects  of  Thickness 

The  nonlinear  effects  of  tliickness  on  unsteady  loads  liave  been  dis¬ 
cussed  in  connection  with  supersonic  strip  tlieory,  and  a  corresponding 
effect  exists  for  slender  bodies  due  to  the  axial  flow  perturbation  interaction 

with  the  crossflow  perturbations.  These  second-order  theories  have  been 

86 

pioneered  by  Lighthill  and  Van  Dyke  for  tlie  steady  state  case.  Van  Dyke 
has  introduced  the  concept  of  a  "hybrid"  theory,  wherein  the  first-order 
crossflow  solution  is  added  to  the  second-order  axial  flow  solution,  and  tlie 


exact  isentropic  pressure  and  surface  Ijoundary  conditions  are  used.  Van  Dyke 
argues  that  approximations  to  the  pressure  coefficient  and  boundary  conditions 


are  unnecessary  even  though  justifiable  on  an  order  of  magnitude  basis.  Tobak 

12s 

and  Welirend  have  extended  this  concept  to  unsteady  rigid  motions  of  pitching 
and  vertical  acceleration  and  present  plots  of  stability  derivatives  for  cones 
including  the  steady'  state  results  and  comparisons  with  exact  results,  and  with 
Newtonian  tlieory  for  large  values  of  Tliis  does  not  constitute  a  second- 

order  crossflow  tlieory;  liowever,  the  use  of  the  exact  pressure  and  surface 
boundary  condition  involves  the  axial  flow  potential  which  is  taken  from  tlie 
second-order  theory  of  Van  Dyke.  The  results  for  a  cone  show  a  reduction 
in  magnitude  of  all  normal  force  coefficients  relative  to  simple  slender -body 
tlieory,  as  increases.  We  now  consider  the  second-order  crossflow 

problem. 

87 

Revell  has  carried  out  a  second-order  unsteady  crossflow  potential 
solution  based  on  slender -body  theory  as  the  starting  point  in  an  iteration 
procedure.  This  solution  refines  both  the  boundary  condition  and  the  pressure 
coefficient  expressions  to  obtain  consistent  approximation  throughout.  The 
disturbance  potential  $  is  split  into  tlie  first-order  slender-body  potential  0  , 
plus  a  second-order  potential  4^  .  Thus,  in  a  body  axis  cylindrical  coordinate 
system  (x,  r,  P)  inclined  at  an  angle  a  the  fluid  motion  is  defined  by  the  dis¬ 
turbance  potential  such  that,  following  Ref.  86, 


(Z04) 
(Z05) 

The  potentials  0  and  '  are  resolved  in  terms  of  the  axial  and  crossflow 
potentials  defined  liy  the  relations 


^  ^  '^1 


(^06a) 


'J'  f  'i  ,  sin0 
0  1 


(Z06b) 
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Now  the  nonlinear  potential  equation  in  fixed  axes  which  is  to  be  satisfied 
is[(See  Ref.  38,  p.  Z,  Eq.  (1.1.9b)J 


a^v‘^0  =  (0/9t  +  yV  n  ■  V  )  (VO)”^  +  Q 

c.  tt 


(Z07a) 


wlier  e 


(7-1) 


n  +  (Mi-) 

T 


(VO)  -  1 


(ZOTl)) 


SiilDstitiition  of  Eqs.  (^04)  to  (Z06)  into  Eq.  {i.07),  expressing  the  time  deriva¬ 
tives  in  the  moving  axis  system  (3/d  t  =  9/3t-Vj^*  V),  .  and  approximating  the 
nonlinear  terms  with  values  computed  from  the  first-order  potential,  0,  leads 
to  the  following  equations 


Or 


+  (1/r)' 


Or 


-  hJK,r) 


(Z0  8a) 


'  irr  +  (1/r)  -(l/r^)Vj  =  Ai(x,r)  .  (208b) 

The  solution  of  Eqs.  (208)  lias  l:)een  carried  out  in  Ref.  87,  wliere  corrections 

to  stability  derivatives  have  been  derived  for  the  low  reduced  frequency  liar- 

nionic  case.  The  results  show  a  reduction  of  magnitude  of  stability  derivatives 

for  cones  as  M^.  increases  which  are  somewhat  greater  than  those  shown  l)y 
b  y  ^ 

’r()l)ak  and  Welirend,  partly  because  tlicy  have  not  computed  the  second-order 
crossflow  potential  contrilnition,  and  partly  due  to  the  rather  different  ap - 
proaclies  involved.  The  si  ender -l^ocly  iteration  j)rocedure  is  considered  valid 
only  for  «<.  0.  7,  beyond  which  an  ai:)parent  reversal  of  trends  is  indicated 

l)y  the  tlieory.  I'lie  second-order  slender-body  solution  descril)ed  alcove  is 
similar  to  tlu'  steady  state  second  order  slender-body  crossflow  solution  devel¬ 
oped  by  Ijighthill;^  th  c  r  (To  rc' ,  the  comparisons  with  more  exact  steady  flow 
theorif's  discussc'd  by  Van  Dvke  are  pertinent.  VVe  now  consider  the  so-callcrl 
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"first-order  tlieory"  (as  opposed  to  "first-order  slender-body  theory"  discussed 
in  th.e  preceding  section). 

The  basic  feature  of  the  first-order  problem  is  that  it  retains  the  wave 

equation  character,  'which  for  the  unsteady  case  in  body  axis  coordinates  may 

be  derived  from  Eq.  (4.  1.  15)  of  Ref.  87,  and  which  agrees  with  the  steady  state 

8  G  *  85 

results  of  Van  Dyke  and  Lighthill.  We  have  for  the  crossflow  potential 


+  (  -  (  l/r 


XX 


xt 


tt 


(209) 


subject  to  tlie  condition 


v>, 

r  r 


-  w  (x,t) 
R  ^ 


and  the  initial  conditions  (for  supersonic  flow) 


0  ^(x,  r ,  t)  =  0 

as  x 0 

0^  (x,  r,  0,  t)  =  0 

X 

and  the  evanescense  at  infinity 


(7  10) 


(211) 


Urn 

r  — >  oo 


0  j(x,  r,  P,  t)  =  0 


(212) 


'Fhe  steady  state  solutions  are  known  and  will  not  be  discussed  (see  Liepmann 

40 

and  I\oshko,  for  example).  T’le  general  unstead  y  probl  em  was  first  solved 

88  89 
by  Dorrance,  althougb  this  solution  w.is  criticizerl  by  Miles  on  the  basis 

(h.it  slender-body  theory  is  the  only  linear  theory  \\^iich  is  consistent  on  the 

bnsis  of  order  of  mrignitude  conditions.  If  is  a  thickness  ])arameter,  the 

coiiditions  are 


for  k  =  0(  1):  8  <.^1  and  I 

for  k»  1:  k8'^<l  and  k^M^o‘^«l 


(213a) 


(213b) 


(See  also  Ref.  38,  Ch.  IZ.  )  Similar  criticism  of  the  steady  flow  analysis  of 
.  90  ...  91- 

Tsien  was  given  by  L,ighthill.  It  is  thus  not  clear  that  slender -body 
theory  can  be  improved  upon  without  going  to  a  second -order  theory.  How¬ 
ever,  it  does  seem  possible  that  the  first-order  theory,  which  includes  tlie 
initial  conditions  on  x  and  improves  accuracy  for  high  frequency  cases,  might 

provide  a  better  starting  point  for  iteration  procedures  as  suggested  by 
87 

Hevell.  Another  unexplored  possibility  might  be  to  \ise  Dorrance's  solution 

in  conjunction  with  a  hybrid  theory  as  proposed  by  Van  Dyke^^  and  by  Tobak 
7  Z 

and  Wehrend,  which  requires  the  axial  flow  perturl^ation  field  to  at  least 
second  order. 

For  the  purposes  of  aeroelastic  analysis,  it  would  be  desirable  to  have  a 
systematic  procedure  for  treating  flexible  lifting  bodies.  The  slender -body 
treatment  given  in  Section  II.  C.  1  is  suggested  as  reasonable  for  range  where 
the  conditioTis  of  Eqs.  (213)  are  satisfied.  For  hy^jersonic  speeds,  as  has  al¬ 
ready  been  noted  in  Section  II.  B.  4,  many  other  considerations  enter.  Recently 
93 

Bond  and  Packard  have  treated  the  first-order  problem  for  a  flexible  body 
for  which  w(x)R  (x)  is  given  by  a  polynomial.  In  their  basic  analysis,  appar¬ 
ently  unaware  of  any  earlier  work,  they  solve  the  identical  problem  treated  by 
88 

Dorrance.  They  show  an  increase  in  stability  derivatives  for  rigid  motion 
with  increasing  Mach  numl)er  for  a  parabolic  nose  l:>ody  of  revolution,  as  com¬ 
pared  to  first-order  slender  body  ("apparent  mass"  theory  in  the  notation  of 

Rc'f.  93).  Tliis  suggestion  is  contrary  to  all  inferences  from  second-order 

89-87  7  Z 

tlieories,  botli  steady  and  unsteady,  the  hyl)rid  theory,  or  the  limiting 

7Z 

Newtonian  values.  The  experimental  data  cited  in  Ref.  93  are  not  convincing, 
in  that  tliere  are  apparent  shifts  in  the  bevels  of  experimental  values  obtained 
in  different  f^icilities  which  cover  the  traiisonic  as  opposed  to  the  supersonic 
r('gime.  It  is  suggested,  thc-refoiu',  that  for  all  supersonic  Macli  numbers, 


-This  can  be  justified  based  o>i  Miles's  extension  of  the  work  of  Warrl 
to  unsteady  flow  (  sec  Ib'f.  38  for  a  bibl  iograpiiy). 
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that  first-order  slender -body  theory  probably  represents  an  upper  bound  on 
stability  on  derivatives  such  as  damping  in  pitch,  etc.  ,  and  that  for  high  Mach 
numbers,  each  section  of  the  body  is  regarded  as  defining  some  range  of 
for  whicli  tlie  appropriate  hypersonic  theory.  This  is  especially  so  in  view  of 
the  importance  of  blunt  nosed  configurations.  For  subsonic  flows,  however,  it 
is  possible  that  the  complete  first-order  theory  miglit  be  an  improvement  over 
slender-body  theory,  and  a  numerical  solution  for  the  doublet  strength  distri¬ 
bution  as  with  lifting  surface  theory  would  be  the  indicated  procedure.  In  this 
connection,  both  Dorrance’s  work  and  that  of  Bond  and  Packard  have  used  the 
slender-body  estimate  of  the  local  source  strength  and  are  not  complete  to  the 
extent  that  the  doublets  in  the  Mach  forecone  of  each  point  on  the  body  are  not 
considered  in  the  evaluation  of  the  doublet  strength  distribution. 

D.  Panel  s 

1.  Two-Dimensional  Panels 

The  two-dimensional  flat  panel  has  been  considered  in  Ref.  94;  we 

review  tliis  derivation  of  the  AICs  in  tlie  present  notation.  From  the  analysis 

95 

of  Nelson  and  Cunningham,  '  the  complex  upper  surface  perturbation  pres¬ 
sure  for  harmonic  motion  can  be  written  in  dimensional  form  as 


p(x)  -  (oV^/?)|li  +  (Zik/  c)jjM^  -  Z)/3^Jh 

X 

(4k^/c^)  j^h  lQx  -  ^)/^de| 


(214) 


where  the  geometry  is  shown  in  Figure  14  for  the  example  of  a  single  span 
panel,  and  the  influence  term  is  given  Ijy 


L(u)  =  cxp(-i'M u)-^  [j +  Z)/ZP  J  Jq(u)u/m) 

+  i(  ZM/8'^)J  u/m)  +  (M^/ZP^)J^((»’u/\4)|- 


(215) 
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Generalizing  Eq.  (214)  to  a  column  matrix  permits  substitution  into  Eq.  (217), 
we  obtain 


[s]  W  -  -(Ax/6k^  (&■}  +  (2ik/c)[(M^  -  2)/B^]{h} 

X  \ 

+  (4k^/c^)|h  l[(x  - 


(218) 


It  remains  to  relate  {_h']’  and  -^^hLd  4  I  to  «  ,  which  can  be  accomplished 
by  finite  difference  methods.  It  is  convenient  to  write 


[S]  =  fhl]  *  f^z]  *  [Sz] 


(219) 


wliere 


[^hl]  ^  -(ix/Sk^)  (k’} 

[S,z]  {■>}  =  (2iAx/ck)(^(2-M^)/p^j|hj. 
fh3]  W  =  -(4Ax/Pcb[(hLcl 


(220) 


(221) 


(222) 


O’hV'ious  I  y 


=  (2i,\x/ck)  ((2-M^)/°-''^'‘ 


(223) 


which,  for  the  single  span  example  panel,  after  clioosing  Ax/c  =  0.  1  (i.  e.  ,  tlie 
oxampde  j^anel  lias  nine  control  points),  Ijecomes  the  following  ninth -order 
constant  imaginary  diagonal  matrix 


79 


(ZZ9) 


ycix  =  (Ax/Z4)  (Vy^  +  i^y^^  -  5y^  f  y^) 


ydx  =  (Ax/Z)  (y^  +  2y^  +  .  .  .  +  Zy^^  ^  -f  y^^) 

MAX//4)  (-y^  + 


(Z30) 


Those  expressions  can  be  combined  to  give  all  of  tlie  required  integrals.  A 
tyj)ical  one  is 


j  li(?)L  -  ?)/c]  (1^  =  Ax  (31/A4)hjL(0.  7)  +  (5/6)  I.^MO.  6) 


+  {A5/A4)h3M0.  5)  +  h^MO.  4) 


+  ha..(0.  3)  +  li,  L(0.  i)  +  (Z5/34)li^l40.  1) 

DO  I 


+  (1/Z)1.  1.(0)  -  (l/Z4)h  M-0.  1)  (Z3I) 

H  9 


Thert’fore,  the  int<.*gr ation  matrix  |c^  for  the  example  panel  with  Ax/c  = 
0.  1  becomes  ’he  following  nintli-order  complex  matrix: 


HI 


Tiie  complete  ninth-order  AlC  matrix  for  tlie  single  span  panel  follows  by  sum¬ 
ming  Eqs.  (^^4),  and  (Z3Z)  according  to  Eq.  (219).  The  extension  of  the 

foregoing  to  multiple  span  two-dimensional  panels  is  straightforward.  The  AICs 
for  a  two  span  panel' having  18  control  points  are  also  shown  in  Ref.  94. 

2.  Th ree -Dimensional  Panels 

The  AICs  for  a  three-dimensional  panel  have  not  been  obtained  spe¬ 
cifically,  The  two-dimensional  analysis  provided  the  simplification  of  a  non¬ 
singular  kernel  in  the  integral  term;  the  three-dimensional  analysis  requires 
the  more  complicated  approaches  of  supersonic  lifting  surface  theory  discussed 
previously  in  Section  11.  B.  3.  In  principle,  these  methods  are  directly  applica¬ 
ble,  but  there  are  practical  numerical  difficulties  because  a  large  number  of 
downwash  control  points  are  needed  to  provide  an  adequate  basis  for  numerical 
solution  for  the  pressure  loading,  in  comparison  witli  the  requirements  for  anal¬ 
ysis  of  wing  modes,  for  example.  This  follows  from  the  fact  that  we  must  have 
a  reasonable  number,  say  six,  control  points  per  half  wave  of  deformation 
(cf.  ,  Ref.  56  in  connection  witli  tlie  box  method).  For  example,  if  we  consider 
the  isolated  panel  of  Fig.  15  to  deform  in  modes  of  the  form  (e.  g.  ,  for  simply 
supported  edges) 


N 

r  s 


h(x,  v)  =  )_  ^  sin(rTTx/a)  sin(  STTy/b) 

r  =  l  s  =  l 


(<>33) 


then  we  require  36  N  N  downwash  control  points  to  obtain  adequate  results. 


r  s 


If  it  is  necess;\ry  to  consider  the  interference  of  N  panels,  then  36  N  N  N 

P  p  r  s 


Ijccomes  tlie  requitu'd  total  number  of  control  points.  By  virtue  of  this  large 
nuiiiber  of  control  pv)inls,  we  thus  conclude  that  the  lifting  surface  theories 
treated  in  Section  11.  IT  3  are  fcasil^le  though  computationally  more  time  con¬ 
suming  (and  expensi'^'  than  most  lifting  surface  studies.  This  application  of 
lifting  surface  tlmorv  to  panel  flutter  investigations,  which  does  not  appear  to 
have  been  made  thus  f.n*,  deserv'cs  further  research.  The  development  of  the 

AICs  follows  the  clevel or^ment  of  Section  II.  B.  3,  once  the  matrix  I  lias 

'  I  nmj 


been  computed,  although  we  note  that  we  must  consider  tlie  upper-  and  lower - 
surface  panels  separately,  and  thus  iiave  only  half  of  the  loading  normally 
developed  for  a  given  downwash  distribution  assuming  that  the  effect  of  tlie  in¬ 
ternal  air  may  be  neglected. 


Fig,  15.  Three -Dimensional  Panel  Geometry 


We  shall  not  consider  the  cylindrical  panel  aerodynamic  problem  in  detail 

except  to  note  as  an  example  that  the  supersonic  problem  has  been  formulated, 

97 

but  not  completely  solved,  by  Holt  and  Strack.  Tlieir  solution  has  been  car¬ 
ried  out  by  Laplace  transform  methods  (transforming  the  chordwise  coordinate) 
for  quasisteady  motion.  Their  treatment  of  the  bounda ry  condition  a t  the  cylin¬ 
drical  surface  is  actually  more  rigorous  than  the  treatments  mentioned  in  Sec¬ 
tion  II.  C.  Z  for  the  slendor-body  problem  for  the  low  frequency  case,  l)ut  be¬ 
cause  they  choose  to  neglect  the  potential  time  derivative  teri7i  in  the  calculation 
of  the  pressure  perturbation,  the  analysis  is  not  an  altogether  consistent  low 
frequency  theory.  This  first-order  analysis  should  be  extended  to  arbitrary 
frequencies,  and  the  estimation  of  second -order  effects  must  be  considered. 

No  attempt  as  yet  has  been  made  to  devcloj:)  the  AlCs  for  cylinders. 
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3. 


Viscous  Effects 


Jf  tiic  longiiuflinnl  and  trai^sverse  wav^c  lengths  of  panel  modes  be¬ 
come  comparable  in  magnitude  to  the  local  boundary  layer  thickness,  then 

'  •  47 

viscous  effects  must  be  considered.  We  have  already  noted  Miles'  work 

in  this  area  in  connection  with  the  discussion  of  vorticity  effects  in  Section 

II.  B.  c.  The  boundary  layer  acts  as  a  nonuniform  free  stream  upon  which 

small  inviscid  disturbances  are  superimposed.  One  of  the  features  of  this 

class  of  problems  is  that  the  pressure  gradient  normal  to  the  surface  is  not 

.  .  77  . 

negligible  as  is  the  case  in  classical  thin  boundary  layer  situations.  This 
assumption  is  the  cornerstone  of  the  justification  for  calculating  pressure  dis¬ 
tributions  from  inviscid  aerodynamic  theory.  The  extension  of  Miles'  ap¬ 
proach  to  consider  the  viscous  disturbances  due  to  surface  motion  for  laminar 
boundary  layers  can  follow  along  the  lines  of  classical  hydrodynamic  stability 

9  8 

theory.  Once  this  extension  is  made,  the  AICs  will  follow  as  before. 
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III.  OSCILLATORY  AICs  FOR  GUSTS 


Recent  advances  in  the  mathennatics  of  random  processes  have  provided 
the  basis  for  statistical  solutions  to  engineering  problems  heretofore  analyzed 
by  rather  arbitrary  deterministic  methods.  Examples  of  interest  to  the  aero- 
elastician  include  gust  stresses,  fatigue  life,  and  servomechanism  performance 
in  a  turbulent  atmospheric  environment.  For  the  analysis  of  these  problems  it 
is  necessary  to  know  both  the  statistical  properties  of  the  forcing  function  (i.e.  , 
power  spectral  density  of  the  turbulence)  and  the  frequency  response  character¬ 
istics  of  the  system  transfer  function  (e.g.  ,  the  stress  at  a  point  due  to  a  har¬ 
monic  gust  with  a  given  wave  length).  The  frequency  response  analysis  for  a 
restrained  system  has  been  briefly  indicated  in  Section  I.  Methods  of  fatigue 
life  analysis  or  servomechanism  noise  analysis  are  beyond  the  scope  of  the  pre  s - 
ent  treatment,  but  having  illustrated  the  need  for  harmonic  gust  analysis,  we 
may  proceed  to  consider  the  derivation  of  oscillatory  AICs  for  a  harmonic  gust. 

We  consider  first  the  case  of  strip  theory.  We  take  the  derivation  from 
Ref.  3.  In  Section  I,  the  gust  AICs  were  defined  by 

If  1  =  pvw  b  s  fc  1  I  w  / w  I 

1  gi  g  1-  [  g]  I  g  g  I 

where  w^  is  the  gust  downwash  velocity  at  each  spanwise  station;  the  column 

Iw  /  W  1  becomes  a  unit  column  matrix  for  a  uniform  (i.e.,  two-dimensional) 
g  gl 

gust  field.  We  have  noted  this  definition  is  completely  general,  being  equally  as 
applicable  to  a  lifting  surface  theory  as  to  a  strip  theory.  However,  in  the  case 
of  a  strip  theory,  the  AICs  take  on  the  simplified  partitioned  form.  For  ex¬ 
ample,  the  gust  AICs  for  a  two  strip  wing  appear  as 


(Z34) 


Eacli  partition  is  of  size  (Z  x  1)  since  two  control  points  are  necessary  on  each 
strip  if  it  has  two  flexible  degrees  of  freedom  (assuming  a  rigid  chord,  and  no 
control  surface). 


86 


In  order  to  derive  tlu'  strip  AICs  we  assunie  that  tiic  gust  lift  and  moment 
about  the  airfoil  quarter  cliord  are  known.  i  lie  equivalence  between  the  given 
loads  and  the  control  point  f{3rces  is  shown  in  h'ig.  16.  We  note  the  sign  con¬ 
vention  is  chosen  anbitrarily,  only  consistency  being  required;  if  the  gust  is 


Fig.  16.  Given  (a)  and  Replacement  (b)  Force 
Systems  and  Geometry 


specified  as  an  upwash,  a  positive  force  acts  upward,  and  if  the  gust  is  specified 
as  a  downwash,  a  positive  force  acts  downward.  From  Fig.  16  we  note  the 
equivalence 


F 


+ 


L 

g 


(235) 


which,  in  matrix  form, 


dF 


2g 


reads 


f 

F, 

L  -1  M  /d 

Ig 

►  -  < 

g  g 

-M  /d 

2g 

g 

(236) 


(2  37) 
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In  order  to  continue  the  derivation  it  is  necessary  to  have  the  relationships 

between  the  lift  and  moment  and  the  gust  velocity*  To  this  end  we  select  the 

99 

incompressible  solution  of  Sears  who  gives  the  lift  and  moment  (when  cor¬ 
rected  for  sweep)  as 


L 

g 


ZircosApVw  bAy $(k)e:cp(-ikx/b) 

g 


(238) 


M  =  0  (239) 

g 

where  the  Sears  function  is  given  in  terms  of  the  Theodor  sen  function  C(k)  as 

$(k)  =  [jptk)  -  ijj(k)|c(k)  +  ij^'k)  .  (240) 


The  first  sweep  correction,  the  factor  cosA,  adjusts  the  two-dimensional  lift 
curve  slope;  the  second  sweep  correction,  the  factor  exp(-ikx/b),  accounts  for 
the  phase  difference  in  the  gust  maximum  amplitude  reaching  the  leading  edge  of 
each  wing  strip.  Substituting  Eqs.  (238)  and  (239)  into  Eq.  (237)  yields  the 
force-gust  velocity  relationship 


Ig 


=  2ncosApVw 


fbAy$  (k)  exp  (-ikx/b)l 

4  °  / 


[241) 


which,  by  comparison  to  Eq.  (3),  yields  the  harmonic  gust  AICs  for  the  j'th 
strip 


2Trcos  A-n 


|(b./b  )(Ay  /  s)  $  (k  )  exp  (-ik  x  /bj) 

J  J  J  J  J 


(242) 


where  the  local  reduced  frequency,  k., 


is  based  on  the  local  semichord  length 


k  =  (I'b./V 

(243a) 

J  J 

=  (2TTV/X)(b./V) 

(243b) 

=  2  n  b^  /  X  , 

(243c) 
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where  X  is  the  gust  wave  length.  As  in  the  case  of  motion,  once  the  AIC 
partitions  are  obtained  for  each  strip,  the  surface  AICs  are  formed  by  assem¬ 
bling  the  partitions  according  to  Eq.  (234). 

N 

We  may  consider  the  formulation  from  lifting  surface  theory  more  briefly, 
since  it  is  only  necessary  to  replace  the  downwash  from  motion  in  Eq.  (55)  with 
the  gust  downwash  matrix  jw  exp(-ik.x./b.)  |  where  the  x.  are  now  the  chord- 
wise  coordinates  of  the  control  points.  Then  the  harmonic  gust  AIC  matrix  from 
lifting  surface  theory  becomes 

l^g!  -  ^  l^iml  [^Ll 

where  the  matrices  and  I  A^  ,  and  the  scale  factor  K  are  defined  in 

nm  J  I  nm  J 

Section  II.  B.  1  for  subsonic  flow  and  in  Section  II.  B.  3  for  supersonic  flow,  and 

the  exponential  matrix  is  a  diagonal  form.  A  similar  treatment  is  possible  in 

the  case  of  slender  body  theory  by  setting  w  =  W  exp(-k  x/b  )  in  Eq.  (  180)  of 

g  r  r 

Section  II. C  and  proceeding  according  to  the  method  of  that  section. 


IV.  TI^ANSIKNT  AlCs 

TIu'  rur rent  status  of  transient  AK/S  permits  only  a  few  preliminary  re¬ 
marks.  For  tlie  monu'iit,  we  shall  cliseuss  somt'  promising  means  availabh*  to 

\ 

utilize  our  oscillatory  solutions  in  deriving  useful  transient  results,  rather  than 
st'eking  transient  solutions  pi'r  sc. 

For  list'  in  dynamic  stability  anti  control  investigations,  a  Taylor  series 
format  may  be  obtained  for  the  transient  AlCs  for  motion.  The  series  may  be 
written 

|f|=  (qs/c)([cj  |h|+KJ|hc/v} 

where  all  matrices  thus  defined  are  real.  In  the  case  of  harmonic  motion,  Eq. 
(Z45)  n^ust  agree  with  Eq.  (Z),  at  least  in  the  frequency  range  of  interest,  as  we 
have  already  shown  in  the  limit  of  zero  frequency  [i.  e,  ,  Eq.  (5)]  .  Hence 

Ksl  +  Kohl  -  Ko-^hl 

^  -  Z^lics/S)  |Cj^(k^.))  (^46) 

where  we  emphasize  the  functional  dependence  of  [  ]  on  the  reduced  fre¬ 

quency.  The  equality  of  matrices  in  Eq,  (Z46)  requires  the  equality  of  each 
element 

(Ss)ij  ^  (SDh)i.i  -  ''r("/^)"(SD^h)ij 

.ikj(c/b,.)^(c,^^3,;)i^  +  •••  =  .kJ(cs/S)  (C,(k_.))... 

(247) 
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Eq.  {1^1)  is  tlie  basis  of  solving  for  Iho  unknown  matrices  of  Eq.  (^4S).  !)(»- 

pending  on  tlie  point  of  truncation  of  the  series  on  the  left-hand  side  of  Eq. 

[141),  a  sufficient  numb(‘r  of  reduced  frequencies  may  be*  chosen  to  yield  a  de¬ 
terminate  system  of' sinuiltaneou s  equations  for  the  ufikiiown  elements.  For 
example,  choosing  only  one  frequency  leads  to  equating  the 

real  and  imaginary  parts  of  both  sides  of  Eq.  (Z47)  (noting  that  is  complex, 
and  assuming  has  been  determined  previously);  choosing  two  frequencies 

leads  to  two  siTiuiltancous  equations  in  C,  and  C,  ^3,  ,  and  two  more  simul- 

taneous  equations  in  C,  and  C,  ^4,  ,  from  whicli  all  four  coefficients  may 

hD  h  hD  h 

be  calculated.  We  anticipate  transient  AICs  derived  in  this  manner  to  have 
practical  value  in  the  study  of  the  dynamic  stability  and  control  of  slowly  ma¬ 
neuvering  vehicles  (i.  e.  ,  tlie  maneuver  may  be  characterized  by  a  low  reduced 
frequency). 

A  more  general  approach  which  is  not  subject  to  the  frequency  restriction 
of  the  above  method  is  based  on  the  convolution  formulation.  Here,  we  may  de¬ 
fine  transient  AICs  by  the  following  expression  [cf,  ,  Ref.  34,  Eq.  (5-370)] 


|f|  =  (qS/c)  lc,,(v(t-g/bj]|h(t„)| 


If  we  assume  tin-  following  exponential  approximation  for  each  element  of  the 
transient  AIC  matrix 


c,  (Vt/b  )..  =  (C,  ).. 

h  r  ij  hs  ij 


L 

k 


b. .  exp(  -  P  • Vt/b  ) 
ijk  ‘  ijk  r 


(^49) 


we  are  in  a  position  to  alig!i  the  harmonic  case  of  Eq.  (Z4H)  with  the  oscillatory 
case  of  Eq.  ( Z) .  The  form  chostni  in  Eq.  (Z40)  leads  to  asymptotic  convergence 
to  the  static  case.  In  orrler  to  evaluali'  Eq.  (Z48)  for  harmonic  motion,  we  let 
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and  st't|^h(l)|  =  lU' j^h  fxp  (i«l)|  .  Nolin^  also  that|^F(t)  J-  =  Rej^F  oxp 

(i(iit)|.,  wt'  find  the  in( ration  loads  to  the  following  result  for  the  complex 
force  amplitude: 


f1-  (qS/c)|c,^(l  - 


(/.50) 


Again  the  conformity  between  Eq.  (^50)  and  Eq.  {Z)  requires  equality  of  each 
element.  Hence 


V 

^  1  -  iR...  /k  “ 

k  ijl< 


,  C  (k  ).. 

=  ZkiZs/S)  p-  pj 
'Ss’ij 


(251) 


As  before,  the  point  of  truncation  of  the  series  determines  the  number  of  re¬ 
duced  frequencies  necessary  to  solve  fur  the  coefficients  of  the  exponential 
series;  each  additional  term  (having  two  coefficients)  requires  one  additional 
frequency  [noting  again  Eq.  (^iSl)  is  complex]  .  For  example,  a  two-term 
series  (perhaps  a  practical  maximum)  requires  two  frequencies  to  solve  for 
the  four  coefficients  in  each  element,  b..,,  B..,,  b..,,  and  R..,.  The  reduced 
frequencies  may  be  cliosen  to  cover  a  wide  range  without  restriction,  although 
good  low  frequency  In-luivior  is  still  desirable.  A  similar  flcvelopmcmt  to  the 
present  one  is  possible  to  obtain  transient  gust  AlCs  from  the  oscillatory  gust 
AICs.  Current  investigations  have  not  progressed  to  this  point,  nor  have  they 
progressed  sufficiently  to  assess  practical  difficulties  in  applying  Eq.  (^i51) 

(e.  g.  ,  wliat  are  the  consequence's  of  an  I’lement  .j  vanishing?).  Investi¬ 

gations  of  transient  AICs  per  se  are  ci’rtainly  not  to  be  discouraged,  but  per¬ 
haps  should  await  demonstration  of  the  utility  of  tlie  foregoing  possibilities. 


Tlie  problem  of  obtaining  the  transient  rc’sponse  of  a  system  should  also 
be  discussed,  since  the  response s  to  discrete  gusts  and/or  blast  waves  areprob- 
lems  of  continuing  concern  to  the  aeroelastician.  We  mention  hc're  an  alterna¬ 
tive  to  tlie  integration  of  the  equations  of  motion  to  obtain  the*  time  history.  We 
recognize  Eq.  (17)  as  a  Fourier  transform  of  the  transient  resfionse  problem. 
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If  the  solution  for  the  response  matrix  |  |  is  carried  out  for  several  values 

of  fr(‘quency,  tiie  transient  response  |  (h^( t)|  may  be  obtained  by  means  of  the 

nunurical  inverse  Fourier  transform  techniques  of  Lanezos  (Ref.  37,  p.  Z63). 
I’lie  calculation  procedure  is  given  by 


h^(t)  = 


h  (k/ZT)  exp  (TTikt/X) 


(^5Z) 


K  =  -  k 

max 


where ‘T  is  a  representative  decay  tiirie,  selected  so  that  jhj^(t)j— 0  whim  t  <-  T 
and  t  >  +  T,  is  taken  as  three  times  the  largest  participating  nat¬ 

ural  frequency  of  the  system.  The  investigations  of  L.  G.  Johnson  and  T.  E. 
Stenton  of  North  Ainerican  Aviation,  Inc.  ,  (private  communication)  have  shown 
that  it  is  not  necessary  to  evaluate  the  AICs  at  each  of  the  frequencies  required 
in  Eq.  (Z5Z)  since  the  AICs  are  smoothly  varying  functions  of  frequency  and 
may  he  readily  interpolated  by  means  of  any  of  several  available  numerical 
interpolation  methods.  Transient  AICs  may  also  be  obtained  directly  by  this 
method  for  use  in  Eq.  (Z48)  since  the  Fourier  transform  technique  is  an  alter¬ 
native  to  the  (’xponential  approximation  of  Eqs.  (Z49)  to  (Z51). 


V.  AN  EMPIRICAL  WEIGHTING  MA'nUX  FOR  USE  WITH  THEORETICAL  AICs 


Historirnlly  s[>('akin^,  both  at* rod yiin m ic  anti  structural  onginocrs  liavt* 
utilizt'd  wind  tunnel  data  in  dt'vising  s emi -en  .pi r  ical  mt'thods  of  aoroelastic 
analysis  for  use  in  dynamic  stability  and  load  studies,  whereas  the  fluttt*r 
engint'er  has  lieen  forced  to  rely  on  theoretical  aerodynamic  data.  Tht*  variation 
in  empliasis  between  experimental  and  theoretical  approaclies  is  justified  to 
some  extent  because  of  tiie  range  of  reduced  frequencit's  appropriate  to  eacli 
problem.  However,  the  semi -empirical  quasi-steady  mt'thods  have  the  distinct 
disadvantage  of  being  attempts  at  deriving  tlie  general  from  the  particular,  but, 
correspondingly,  the  theoretical  methods  have  the  disadvantage  of  not  approach¬ 
ing  the  experimental  values  as  the  reduced  frequency  is  made  smaller.  It  is 
sugg('sted  that  a  unification  is  possible  if  the  tlieoretical  data  are  adjusted  to 
match  the  available  experimental  results,  and  thus  the  generality  of  the  original 
formulation  may  be  retained. 

Following  Refs.  100  and  101,  we  propose  thal  a  premultiplying  diagonal 
vvH'ighting  matrix  be  introduced  into  Eqs.  (1)  and  {,1)  to  adjust  the  tlieoretical 
AICs  to  yield  agreemc'iit  with  ex{)«* r imental  data.  The  corrected  static  force 
d  i  s  t  r  i  hu  t  ion  be  comes 

jl'j  =■  (qS/c)  [vV  ]  [C|^J  |h|  .  (Z53) 

and  tht*  corrected  oscillatory  force  distribution  for  motion  l:)econics 

|f|  s|w]|c,J  |h|  .  (Z5.1) 

I'qs.  (^‘^3)  and  {Z54)  indicate  the  form  in  which  tlic  weighting  matrix  is  to  ])e 
utilized  in  aeroelastic  analysis.  If  tlu*  experimental  data  correspond  to  a  rigifl 
untwisted  configuration,  we  set  jh  |  =C2  jxjin  F'q.  (Z53)  and  pose  the*  problem  as 

tlu'  solution  of 

j  pj  z-  qsa  [w]  [('i  J  jx/:|  (ZS'3) 

fo)'  I  \V  j  pivt’n  all  tbr  rcniaintiig  tci'ms  in  fhc  equation.  We  consicU'r  two  pos- 
sibililies  regarding  the  extemt  of  our  (' xpe  r  inumtal  knowledge: 

I.  Only  a  meagf'r  amount  of  rlata  are  available,  e.  g.  ,  the*  lift  curve 
sloi)e,  and  the*  spanwise  and  ehoialwise  location  of  the  aerodynamic  cauiler,  or 


A  vv('aUh  oi  data  ari*  availal)lc,  t*.  ,  IIk'  spaiiwisi*  variati(.)ns  of  lift 

distribution  aiid  (.bore!  wisr  location  of  local  a(' r  od  ynaniic  cc’nicr.  C.isc  Z  pre¬ 
sents  no  difficnlty  sinci*  a  snfficic'iit  nund)cr  of  (.‘qnations  can  bo  written  from 
Fiq.  (ZSS)  to  (h'tt'rmune  all  clcMiHMits  of  |wj.  (Note:  In  this  ease,  the  weighting 
matrix  nccfl  not  bo  diagonal,  and,  dopemding  on  tbc  aerodynamic  theory  used, 
othi'r  formats  may  bo  profo  rahlo.  )  In  Case  1  thi'rc'  are  considerably  more  un¬ 
knowns  to  1)0  found  than  equations,  and,  lu'uco,  wo  again  turn  to  the'  motlioci  of 
least  squares,  althougli  wt*  no*  e  W(*  arc'  hi're  dc'aling  witli  an  unde  rdc'tc'r  mined 
system  in  contrast  to  the*  ovc'rd  etc*  rminc'd  systems  discussed  previously.  Wc' 
propose  our  additional  equations  from  the  rc^quirernent  that  tlie  changes  in  the 
theoretical  load  d istr il^ution  shall  be  as  uniform  as  possible  or,  in  least  squares 
terminology,  thc'  wc'ighted  sum  of  the  squares  of  the  deviations  shall  be  a  mini¬ 
mum,  where  tlie  deviation  is  defined  iis  tlie  difference  of  final  and  original 
weigliting  factors 

j  e  j  =  I  W  -  I  ]  .  (ZS6) 

The  weighting  func:tion  we  choose  is  the  theoretical  load  distril)ution.  (Tliis 
choice  leads  to  equal  weighting  elements  for  the  case  wlicre  only  the  lift  curve 
slope  is  known.  )  The  weighted  least  squares  condition  tlien  becomes 

E  c  ^  =  j  e  I  j  j  e  I  =  <1  minimum  (Z57) 

wlu'i'c  I  lu*  middle'  matrix  factor  is  a  diagonal  form  of, 

=  qsa  |c^J|x/7  I  . 

It  is  now  conveni('nt  to  make  a  distinction  betwec'U  tbc*  weighting  factoi's 
associated  with  the*  numbc'r  of  givc'n  conditions  and  those'  associated  with  the  re¬ 
in  lining  Ic'ast  sfjuares  conclitions.  In  analogy  with  tlu*  prohlc'in  of  redundant 
struetuiu's  wc  lU'fc'r  to  tlu*  given  conditions  as  basic-  c'quations  and  tlu*  corres¬ 
ponding  weighting  factors  as  l)asie  weighting  factors  »  and  tlie  arlditional 
equ.it  ions  .IS  redund.ant  c'quations  with  W  .as  the  corresponding  redundant  wi'igld 
ing  f.utors.  Introducing  this  distinc'tion  into  k'.qs.  (ZS6)  and  (Zb7)  we  h.avc' 


and 


(Z59a) 

(Z59b) 


hb|  M^rl  '  l^r’l  |e,|  =  ^‘  minin'-'"'  •  (‘^^’0) 


If  we  differentiate  Eq.  (Z60)  with  respect  to  each  of  the  redundant  weighting  fac¬ 
tors,  we  obtain  the  redundant  ('qiiations 


(261) 


where,  from  Eqs.  (259) 


’  9^b ' 

aw 

r 

aWj. 

(262a) 


/ 


(262b) 


since  the  eind  are  considered  the  dependent  and  independent  variables, 
respectively. 


We  next  consider  the  basic  equations.  We  may  write  these  in  tlie  form 


or  by  pa  rt  it  ion ing 


(263) 


(26^1) 


where  Q  denotes  a  generalized  force  (lift,  pitcliing  moment,  rolling  monumt, 
ite.)  and  |  ^  j  is  the  given  set  of  expc' r  ime  ntal  conditions  to  be  matched. 

!"»)r  c'xamplc',  if  th(’  tliiU'e  conditions  givem  are  the  lift  curve  slope*  and  aerorly- 
namic  cent<.*r  coordinates,,  then 

(26  5) 
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s  i  mil  1 1  ;i  noon  s  solution  o[  Eqs.  (261)  and  (2(»4)  is  straightforward, 
h'roio  L^q.  (26d) 


(DI-l 


from  which,  by  diffi' l•('ntial  ion. 


Then,  combining  Eqs.  (259),  (262),  and  (266)  into  Eq.  (261)  permits  solution 
for  the  redundant  weighting  factors 


where 


Iw^  =  [a]  -Mb 


= 

I  ^  i  aii;  ^  ^  J  'AMT 


(26H) 


(269a 


B  = 


+  F' 


l  b 


(t)l-l|^(e) 


(269b) 


and  where  the  derivative  matrix  is  given  in  Eq.  (267).  The  basic  weighting 
factors  are  found  by  substituting  Eq.  (268)  back  into  Eq.  (266),  and  the  weight¬ 
ing  matrix  is  finally  found  by  properly  assembling  the  basic  and  redundant 
weighting  factors  into  the  diagonal  format. 

It  has  been  sllown  theit  a  weighting  matrix  can  be  defined  which  provides  a 
means  of  unifying  the  use  of  theoretical  AICs  and  experimental  data.  The  result 
of  tliis  unification  is  that  the  resulting  ompirically  modified  matrices  of  static 
or  oscillatory  AICs  retain  the  generality  of  the  original  theoretical  basis  while 
including  the  experimental  results  as  limiting  values. 


VI.  CONCLUDING  REMARKS 


We  conclude  that  the  status  of  AICs  is  reasonably  well  advanced,  particu¬ 
larly  in  the  case  of  oscillatory  motion  of  lifting  surfaces.  A  multiplicity  of 
theoretical  solutions  for  AICs  has  l)t*en  established,  and  a  measure  of  unifica¬ 
tion  has  been  achieved  for  a  diversity  of  configurations  and  throughout  the  range 
of  Mach  numl)er  of  current  interest;  the  empirical  modification  of  the  theoretical 
AlC  s  suggestc'd  furtlier  extends  the  unification.  I'hc  generality  of  AICs  makes  it 
liighly  desirable  tliat  subsequent  investigations  be  carried  through  to  this  point, 
particularly  for  computational  efficiency.  It  is  hoped  that  the  format  presenJed 
here  can  provide  a  canonical  form  for  continuing  developments. 

A  number  of  recommendations  can  be  made  for  further  study  beyond  those 
already  suggested  in  the  text  .  These  range  from  research  on  fundamental  prob¬ 
lems  to  numerical  implementation  of  established  solutions  in  the  A IC  format .  Of 

particular  impo  rtance  are  developments  in  the  transonic  and  hypersonic  regimes. 

39 

In  the  transonic  connection,  the  work  of  Landahl  has  provided  a  substantial 
basis  for  improved  transonic  calculations.  In  the  hypersonic  connection,  the 
aeroelastician  must  take  his  lead  from  the  field  of  ae rothermodynamics  for  fur- 
tlier  developments.  In  addition,  there  is  an  urgent  need  for  further  study  of 
slender  bodies  and  the  interference  of  bodies  and  low-aspect  ratio  fins.  It  is 
also  recommended  that  a  number  of  harmonic  gust  AIC  solutions  be  established 
based  on  the  lifting  theories  discussed  here.  Einally,  the  possibilities  suggested 
for  transient  aerodynamic  calculations  warrant  further  investigation. 

It  is  difficidt  to  assess  tlie  need  for  expe  rimental  correlation.  Most  of  the 
theories  discussed  here  have  already  obtained  some  measure  of  correlation,  and 
this  correlation  also  applies  when  the  theory  is  cast  into  the  AIC  format.  How¬ 
ever,  one  candidate  for  a  correlation  study  is  the  empirical  weighting  matrix 
discussed  in  Section  V.  Tlic  correlation  in  the  static  case  will  have  to  be 
carried  out  from  force  measurements  on  twisted  models;  the  correlation  in  the 
oscillatory  case  may  be  made  from  flutter  data.  Such  studies  should  show  not 
only  the  accuracy  of  this  method  bat  also  may  indicate  means  of  finding 
a  more  gcmeral  format  for  tlie  weigliting  mat J'ix  than  the  real  diagonal  form 
presently  consider  ed . 
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